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Abstract 


m 


The  behavior  of  soil  subject  to  a  wide  range  of  stresses  is  of  interest  in  the 
design  of  Air  Force  strategic  and  technical  facilities.  For  an  adequate  facil¬ 
ity  design,  it  is  important  to  determine  accurately  the  soil  properties  for 
the  loading  and  environmental  conditions  under  consideration.  As  a  result,  the 
survivability  and  vulnerability  of  AF  facilities  can  be  better  predicted. 

The  objective  of  the  work  reported  here  was  to  explain  the  mechanical  proper¬ 
ties  of  fine-grained  soils  in  terms  of  the  inter-particle  and  inter-aggregate 
effects  which  control  them.  Through  this  study,  structures  of  natural  soils 
can  be  better  described  and  catalogued. 

The  approach  was  to  use  electron  and  polarizing  optical  microscopy  to  study 
changes  in  clay  microstructure  under  different  stress  conditions  by  using  dif¬ 
ferent  stress  paths  in  triaxial  tests. 

New  methods  of  image  analysis  were  developed  for  both  types  of  microscopy.  Use 
of  these  new  techniques  has  shown  that  the  raicromechanics  of  deformation  are 
more  complicated  than  might  have  been  imagined.  Initially,  all  samples  had 
horizontal  preferred  orientation.  In  both  drained  and  undrained  normally- 
consolidated  kaolin:  from  0-8%  strain,  the  specimens  appeared  to  bed  down  with 
an  increase  in  strength  of  anisotropy;  from  8-16%  strain,  this  tendency 
appeared  to  reverse:  and.  at  about  16%  strain,  failure  planes  developed  non- 
conformably.  In  over-consolidated  kaolin:  from  0-8%  strain,  there  appeared  to 
be  very  slight  bedding  down,  followed  immediately  by  the  non-conformable 
development  of  failure  planes. 


1  Glasgow  University,  Glasgow  G12  800.  Scotland. 
University  of  East  Anglia.  Norwich  N  R  4  7  T  J  ,  England. 
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Chapter  1 


I  NTRODUCT ION 


1 .  1  INTRODUCTION 

The  principal  objective  of  this  investigation  has  been  to  observe  the  changes 
in  microstructure  of  clay  soils  as  they  were  deformed:  in  particular,  the 
changes  during  triaxial  shear  in  samples  of  kaolin  have  been  studied  under 
three  conditions:  normally  consolidated  drained  shear;  over-consolidated 
drained  shear;  and  normally  consolidated  undrained  shear.  These  changes  were 
found  to  be  more  complex  than  had  originally  been  anticipated;  and  the  results 
obtained  provide  a  test  set  against  which  micro-mechanical  computations  may  be 
evaluated.  These  results  are  reported  in  Chapters  5,  6,  and  ",  respectively. 

The  changes  of  the  microstructure  during  deformation  were  measured  using  new 
methods  of  image  analysis  for  both  electron  and  optical  microscopy  which  were 
specially  developed  for  this  work.  Both  sets  of  methods  represent  considerable 
advances  over  what  was  available  before:  and  considerable  interest  has  been 
shown  in  this  side  of  the  work.  These  new  methods  are  reported  in  Chapters  3 
and  4. 

Almost  all  of  the  scanning  electron  microscopy  used  the  back-scattered  mode  of 
observation  on  samples  which  had  been  impregnated  and  ground  to  reveal  plane 
cross-sections,  the  work  representing  one  of  the  first  major  studies  using 
this  technique. 

The  mechanical  testing  procedures  are  reported  in  Chapter  2. 

Publications  resulting  from  this  project  are  listed  in  Appendix- 1.  Other 
references  are  in  Appendix  2.  Figures  and  Tables  are  placed  at  the  end  of  each 
chapter . 


1.2  SCHEDULE 

During  the  first  year  of  the  research,  the  main  emphasis  was  directed  to  con¬ 
structing  three  medium  pressure  filter  presses  to  a  new  design,  instrumenting 
these  and  the  triaxial  rigs,  tuning  the  specimen  preparation  techniques, 
investigating  and  standardizing  the  optimum  conditions  of  observation  in  the 
scanning  electron  microscope,  and  setting  up  the  systems  to  acquire,  archive 
and  analyze  images. 


The  normally  consolidated  drained  series  of  tests  was  started  first;  and, 
apart  from  a  short  pause  whilst  the  medium  pressure  filter  presses  were 


stripped  and  checked,  the  over-consoi idateu  drained  series  followed 
immediately.  Experimental  work  on  this  series  continued  until  the  final  days 
of  the  project.  The  normally  consolidated  undrained  series  of  tests  was  con¬ 
ducted  in  parallel  using  other  apparatus. 

The  hardware  for  image  archiving  and  analysis  was  expanded  steadily  throughout 
the  life  of  the  project:  and  terming  and  tuning  the  new  algorithms  and  imple¬ 
menting  further  new  ideas  was  also  a  continuous  task. 


1  . 3  HYPOTHESES 

The  materials  studied  were  kaolin;  and  almost  all  the  particles  were  small 
fiat  plates.  It  was  expected  from  earlier  work  (e.g.  Smart  1966b,  Smart  and 
Tovey  1981)  that  most  or  all  of  these  plates  would  come  together  face-to-face 
in  sub-parallel  arrangements  to  form  groups  called  domains.  It  was  also 
expected  ( loc.  cit.)  that,  when  failure  planes  formed,  the  plates  within  the 
planes  would  be  rearranged  to  lie  parallel  to  the  planes.  The  immediate  inter¬ 
est,  therefore,  was  how  the  plates  and  domains  would  behave  during  pre-peak 
shear  deformation.  The  initial  hypotheses  may  be  summarized  as: 

1.  Domains  grow  in  size  and  gradually  turn  into  a  preferred  direction  parallel 
to  the  failure  plane.  Fig.  1.1. a;  this  hypothesis  dates  back  to  an  observation 
of  Smart  (1966c). 

2.  Domains  act  as  rigid  bodies  and  slip  into  closer  packing  without  sig¬ 
nificant  change  in  size  or  internal  packing,  Fig.  1.1. b;  this  hypothesis 
springs  from  McConnachie’s  (1974)  study  of  consolidation. 

3.  Domains  retain  their  general  identity,  orientation,  and  relative  position, 
whilst  strain  is  accommodated  by  minor  changes  in  their  size  accompanied  by 
small  changes  in  particle  packing  as  the  latter  slip  over  one  another,  Fig. 

1. l.b. 

4.  Domains  break  into  smaller  domains;  and  some  of  the  fragments  rotate  into 
new  directions,  Fig.  1.1. c. 

It  had  also  been  suggested  (Smart  1985)  that  deformation  increases  the 
anisotropy  of  isotropic  clays  and  decreases  the  anisotropy  of  strongly 
anisotroipc  clays,  i.e.  that  there  may  be  a  critical  anisotropic  state  just  as 
there  is  a  critical  void  ratio.  The  basis  for  this  suggestion  is  summarized  in 
Fig.  1.2,  which  shows  the  change  in  anisotropy  during  shear  deformation  versus 
the  initial  anisotropy,  the  measurements  were  made  by  polarizing  micro- 
photometry,  c.f.  Chapter  4.  The  figure  also  shows  a  regression  line,  which  is 
significant  at  the  7  %  level.  Since,  the  measurements  were  made  by  various 
workers  at  various  times  using  various  apparatus,  and  since  various  clays  were 
tested  under  various  conditions  in  various  apparatus;  this  level  of  sig¬ 
nificance  is  reasonably  high.  The  implication  is  that  deformation  may  be 
accommodated  by  two  different  mechanisms  in  strongly  and  weakly  anisotropic 
clays,  respectively.  (The  figure  is  based  on  data  of  Andrawes  et  al .  1975, 
Bondarik  et  al.  1970,  McKyes  and  Yong  1971,  Omer  1978,  Shibakova  1970,  Smart 
1966a,  Smart  and  Dickson  1979.) 
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.4  ADMINISTRATION 

Ihe  work  was  funded  by  AFOSR  for  four  ."ears  and.  including  no-cost  extensions, 
the  period  of  study  was  from  1  Sep  until  29  Feb  1992. 

The  work  was  split  between  the  University  of  Glasgow  (Glasgow)  and  the 
University  of  East  Anglia  ( (JEA ) .  The  work  in  Glasgow  then  attracted  additional 
funding,  and  all  the  work  there  became  completely  integrated  as  indicated 
below. 

The  work  had  been  planned  by  Drs. Smart  and  Tovey,  who  were  responsible  for 
supervising  the  work  in  Glasgow  and  UEA,  respectively.  The  following 
Assistants  were  directly  employed: 


Dr. Mark  W. Hounslow 

UEA 

RA.  1  .A 

1 

Feb 

88-31 

Aug 

91 

Mrs. (now  Dr.)  Xiaoling 

Leng 

Glasgow 

RA. 1 .A 

1 

Sep 

88-31 

Dec 

90 

Dr.Olumuyiwa  A.Awoleye 

Glasgow 

RA. 1 .A 

1 

Jan 

91-30 

Apr 

91 

Miss  J.P.Desty 

UEA 

Casual 

1 

Apr 

91-30 

Jun 

91 

Glasgow* 

RA. 1 .B 

1 

Jul 

91-31 

Oct 

91 

Mrs. (now  Dr.)  Xiaohong 

Bai 

Glasgow 

RA. l.A 

1 

Oct 

91-31 

Dec 

91 

Miss  M.D. Martinez 

Glasgow* 

RA. 1 .B 

1 

Nov 

91-29 

Feb 

92 

*  wo  rkeil  at  UEA. 


Mrs.Leng  continued  to  contribute  to  the  project  whilst  supported  by  SERC  from 
1  Jan  1991  until  31  Dec  1992.  Mrs.Bai  was  supported  for  three  years  by  Taiyuan 
University  of  Technology.  As  an  undergraduate,  Miss  Martinez  had  completed  a 
study  of  Intensity  Gradient  Formulae  before  taking  up  employment. 

A  study  of  Hough  transforms  was  undertaken  for  us  by  Mr. (now  Prof. Dr.)  Luciano 
da  F. Costa  and  Dr. Mark  B. Sandler,  King’s  College,  London,  supported  by  FAPESP, 
Brazil;  and  further  image  analysis  was  undertaken  by  Mr.Daisheng  Luo  and 
Dr. James  E.S. MacLeod,  Department  of  Electrical  Engineering,  Glasgow,  supported 
by  the  University  of  Sichuan. 

Babcock  Power  pic,  who  had  already  checked  the  design  calculations  for  the 
medium  pressure  filter  presses,  also  gave  some  specialized  help  in  their  fab¬ 
rication.  SERC  Transputer  Initative  provided  additional  image  analysis  facili¬ 
ties  in  Glasgow;  to  which  RSRE  contributed  some  software. 

Interim  Reports  were  submitted  on  26  Oct  1988,  26  Oct  1989,  8  Jul  1990  (Con¬ 
tinuation  Proposal),  and  28  Jan  1991.  Drs. Smart  and  Tovey  attended  AFOSR  Con¬ 
tractors’  Meeting,  Boston,  1987.  Dr. Smart  made  presentations  to  AFOSR, 
Washington,  and  AF  Geotechnical  Laboratory,  Albuquerque,  1990:  Dr. Tovey  made  a 
presentation  to  AF  Geotechnical  Laboratory,  Albuquerque,  1990.  Dr. Smart 
attended  AFOSR  Contractors’  Meeting,  Albuquerque,  1992.  Dr.Leng  attended  AFOSR 
Contractors'  Meeting,  Boston,  1993. 

As  a  result  of  the  work  funded  by  AFOSR,  Drs. Smart,  Tovey,  and  Hounslow  have 
received  subsequent  funding  from  either  SERC  or  NERC;  Drs. da  Costa  and  Bai 
expect  to  receive  support  from  their  national  bodies;  and  Dr.Leng  has  put  for¬ 
ward  equally  important  proposals. 
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Fig. 


Fig.  1 


1.1  Hypothetical  modes  of  deformation;  two-dimensional  idealization. 
Left:  Original  structure.  Right:  At  failure. 

a.  Top:  Hypothesis  1.  Domains  grow  larger  and  re-align. 

b.  Center:  Hypothesis  2.  Domains  slip  into  closer  packing,  or 

Hypothesis  3.  Particles  slip  within  domains. 

c.  Bottom:  Hypothesis  4.  Domains  break  into  smaller  domains. 

AA 


2  Change  of  Anisotropy  vs.  Initial  Anisotropy  during  shear  deformation 

of  clay  (from  Smart  1993). 


2.1  INTRODUCTION 


This  Chapter  reports  the  preparation  ana  procedures  used  for  the  three  series 
of  triaxial  tests  for  which  microstructurai  measurements  were  obtained: 
normally  consolidated  drained  shear:  over  consolidated  drained  shear;  and 
normally  consolidated  undrained  shear. 

Each  series  of  microstructurai  measurements  required  replicate  samples  at  var¬ 
ious  stages  of  deformation:  and,  in  order  to  produce  series  of  samples  which 
were  homogeneous  (i.e.  having  the  same  internal  microstructure  and  properties) 
and  repeatable  (i.e.  all  of  each  series  the  same),  clay  was  mixed  to  a  slurry, 
consolidated  uniaxial lv  in  specially  designed  filter  presses,  reduced  to 
standard  size,  and  consolidated  triaxial ly.  The  details  of  this  mechanical 
preparation  and  of  the  triaxial  testing  itself  are  given  below.  The  stress- 
strain  curves  for  the  three  main  series  of  tests  will  be  shown  with  the  other 
main  results  in  Chapters  5,  6,  and  7. 


2.2  MATERIALS 

The  clay  used  in  all  the  tests  was  Speswhite  Kaolin  obtained  from  English 
Clays  Lovering  and  Pochin  Ltd.  The  material  had  been  hydraulically  mined  from 
weathered  granite  near  St.Austell,  Cornwall,  after  which  it  had  been  screened 
and  classified  using  various  physical  and  chemical  treatments  which  did  not 
include  grinding.  One  batch  was  used  for  the  drained  tests,  and  a  second  and 
nominally  identical  batch  was  used  for  the  undrained  tests,  see  Table  2.1.  X- 
ray  diffraction  showed  the  material  to  be  predominately  kaolinite  with  Je«s 
than  5  %  quartz  and  illite.  Some  peaks  were  diffuse,  suggesting  the  kaolinite 
to  be  slightly  disordered.  X-ray  mapping  indicated  that  the  potassium  and  iron 
impurities  were  mainly  associated  with  the  larger  particles.  Electron  micro¬ 
scopy  confirmed  the  platy  nature  of  the  particles. 


2.3  FILTER  PRESSES 
2.3.1  Introduct ion 

This  Section  describes  the  medium  and  low  pressure  filter  presses  which  were 
used  for  the  drained  and  undrained  samples,  respectively;  details  of  their  use 
are  in  the  next  Section. 
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2.3.2  Medium  pressure  filter  presses 

The  drained  samples  were  prepared  in  a  set  of  three  medium  pressure  filter 
presses  illustrated  in  Fig.  2.1.  Saturated  clay  is  contained  in  the  upper  part 
of  a  double  acting  piston-cylinder  arrangement,  and  water  (or  oil)  under  pres¬ 
sure  is  supplied  under  the  piston.  Porous  stones  above  and  below  the  clay 
permi .  the  pore  water  to  escape  up  out  of  the  top  cap  and  down  through  the 
piston  rod,  which  is  hollow.  Flexible  tubing  conducts  this  water  to  a  suitable 
receptacle,  which  must  be  above  the  level  of  the  top  cap  to  prevent  the  hollow 
piston  rod  from  drying  out.  U-rings  around  the  porous  stones  prevent  the 
slurry  from  escaping  past  the  stones. 

The  general  procedure  was  as  follows.  First,  the  top  cap  was  removed,  the 
lower  part  of  the  cylinder  was  deaired  through  an  air  vent  plug  in  the  piston, 
the  piston  was  lowered  as  far  as  possible,  the  lower  porous  stone  fitted,  the 
sides  of  the  cylinder  greased,  a  filter  paper  placed  over  the  lower  porous 
stone,  and  the  cylinder  carefully  filled  with  slurry.  Then,  the  slurry  was 
leveled  using  a  special  scraper,  the  top  filter  paper  and  porous  stone  fitted, 
an  insert  laid  over  these,  the  top  cap  replaced,  and  pressure  admitted,  see 
below.  When  the  cake  was  ready,  the  inlet  valve  was  shut,  the  top  cap 
removed,  and  the  cake  extruded  by  careful  use  of  the  inlet  valve. 

These  filter  presses  were  designed  to  produce  filter  cakes  4  inches  (100  mm! 
high  and  5  inches  (125  imn)  diameter  from  slurries  which  were  not  unduly 
aggressive  and  which  were  free  from  strongly  abrasive  constituents:  there  was 
provision  to  add  a  protective  layer  0.02  inches  thick  within  the  bo,  e  should 
this  have  become  necessary.  The  material  is  Type  316  stainless  steel  (see  BS 
9'70 :  lO^O). 

The  base  plate  contains  a  specially  designed  non-return  valve  fitted  at  the 
high  pressure  inlet,  a  safety  valve,  and  a  port  for  a  pressure  transducer. 
Movement  of  the  piston  can  be  monitored  from  the  piston  rod. 

The  stress  calculations  had  been  made  in  accordance  with  BS  5500:  1985  as 
Category  3  constructions,  i.e.  the  welds  would  be  accepted  on  the  basis  of 
visual  inspection,  taking  the  operating  pressure  as  250-psi  (1725  kPa)  for  an 
air- water  pressure  source  and  500  psi  (3450  kPa)  for  pure  hydraulic  pressure. 
The  corresponding  standard  test  pressure  would  then  have  been  1000  psi. 
However,  both  the  original  source  of  material  and  the  original  arrangements 
for  welding  had  to  be  changed  at  the  last  moment,  so  it  was  decided  to  test  at 
600  psi.  The  maximum  operating  pressure,  therefore,  was  150  psi  (1035  kPa), 
which  was  just  above  the  maximum  pressure  available.  It  is  thought  that,  if  a 
high  pressure  source  were  available,  and  if  the  presses  could  be  run  in  an 
enclosed  environment,  then,  by  increasing  the  dimensions,  applying  full 
inspection  and  improving  the  top  seal,  a  working  pressure  of  at  least  1000  psi 
could  be  obtained  with  this  design  of  filter  press. 

In  making  the  stress  calculations,  some  difficulty  was  experienced  in  deciding 
the  proper  uplift  pressure  on  the  top  cap.  The  usual  method  of  calculation 
assumes  Leakage  from  the  bore  radially  outwards  along  the  gasket,  taking  full 
internal  pressure  at  the  bore,  and  assuming  that  this  decreases  uniformly 
until  it  is  completely  relieved  at  the  pitch  radius  of  the  bolts.  This  assump¬ 
tion  is  not  quite  safe  for  a  slurry,  which  might  be  extruded  slowly  enough  for 
it  to  dry  out  and  harden  at  or  beyond  the  pitch  radius,  so  that  fuj i  pressure 
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could  build  up  to  a  greater  extent  than  is  generally  assumed.  Consequently, 
the  underside  of  the  top  cap  was  relieved  with  an  upward-sloping  cut  in  an 
attempt  to  control  the  behavior  of  any  slurry  which  was  extruded.  No  extrusion 
of  slurry  was  reported  during  operation;  but  there  was  leakage  during  the 
proving  tests,  and  these  were  completed  by  tightening  the  bolts  to  a  greater 
torque  than  had  been  calculated.  It  is  not  knov.r.  whether  the  method  of  cal¬ 
culation  was  faulty  (although  it  was  taken  from  BS  5500:  1985  and  had  been 
checked  independently),  or  whether  the  material  properties  of  the  gasket  had 
been  wrongly  estimated,  or  whether  the  relief  had  permitted  the  top  cap  to 
flex  and  rise. 

The  operating  temperature  was  specified  as  7  to  50  degrees  Centigrade. 


2.3.3  Low  pressure  filter  presses 

The  undrained  samples  were  prepared  in  a  set  of  three  low  pressure  filter 
presses  of  the  same  basic  design  as  the  medium  pressure  filter  presses,  but 
made  from  quarter-inch  thick  plastic.  These  presses  had  been  used  earlier  by 
Dickson  and  Smart  (1979),  and  no  trouble  was  experienced  with  them  apart  from 
occasional  minor  breakages  due  to  the  age  of  the  plastic. 


2.4  UNIAXIAL  CONSOLIDATION 

For  the  drained  samples,  approximately  2  kg  of  kaolin  powder  was  mixed  with 
2.8  litres  of  distilled  water  (140  %  moisture  content,  3.0  Liquidity  Index). 
The  slurry  was  whisked  until  smooth  in  a  large  food  mixer,  de-aired  by 
vigorous  shaking  under  low  vacuum,  and  carefully  transferred  to  the  medium 
pressure  filter  presses.  Starting  from  a  small  value,  the  pressure  was  doubled 
daily  up  to  700  kPa.  During  the  uniaxial  consolidation,  the  samples  decreased 
in  height  from  320  mm  to  140  mm,  approximately,  (56  %  strain). 

To  make  one  undrained  sample,  2.5  kg  of  kaolin  powder  was  mixed  with  3  kg  of 
de-aired  water  (120  %  moisture  content.  2.56  liquidity  index).  The  slurry  was 
mixed  in  a  Z-bladed  mixer  for  2  hours  under  a  vacuum  of  approximately  500  mm 
Hg.  The  slurry  was  then  consolidated  uniaxial ly  in  the  low  pressure  filter 
presses,  see  Table  2.2.  The  friction  in  each  press  had  been  estimated  by 
noting  the  pressure  required  to  move  the  piston  when  the  press  was  empty,  and 
the  pressures  applied  were  increased  accordingly. 

In  addition  to  the  main  series  of  tests,  an  additional  series  of  8  normally 
consolidated  drained  samples  was  made  in  which  everything  was  done  in  exactly 
the  same  way  as  in  the  main  series,  except  that  a  gauze  was  drawn  up  through 
the  slurry  after  initial  placement  in  the  filter  press  with  a  view  to  improv¬ 
ing  the  homogeneity  of  the  samples.  However,  the  microstructural  examination 
showed  that  the  wakes  of  the  wires  persisted  in  the  samples  throughout 
uniaxial  consolidation,  triaxial  consolidation,  and  shear  deformation.  It  is 
understood  that  a  similar  result  had  been  obtained  elsewhere,  but  details  are 
lacking.  Towards  the  end  of  the  study,  one  of  the  filter  presses  was  extended 
and  filled  by  sedimentation:  this  appeared  to  be  a  more  promising  method,  but 
it  is  slower,  so  more  presses  would  be  required. 
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2.5  TRIAXIAL  CONSOL I DATIQN 

One  standard  triaxial  test  sample  3  inches  high  x  1.5  inches  diameter  (75  mm  x 
37  mm)  was  cut  from  each  filter  cake.  Although  6  samples  were  cut  at  inclina¬ 
tions  between  30°  and  90°;  all  of  the  main  series  of  samples  were  vertically. 

Standard  triaxial  testing  apparatus  was  used. 

For  the  normally  consolidated  drained  samples,  drainage  was  from  the  top,  with 
pore  pressure  measurement  at  the  bottom.  No  back  pressure  was  used,  because 
obtaining  the  highest  possible  effective  pressure  was  more  important  than 
achieving  perfect  saturation.  The  cell  pressure  was  increased  by  stages  to  700 
kPa. 

The  procedure  for  the  over-consolidated  drained  samples  followed  that  for  the 
normally  consolidated  samples,  except  that  the  cell  pressure  was  finally 
reduced  from  700  kPa  to  100  kPa. 

For  the  undrained  samples,  drainage  was  from  the  bottom,  with  a  plain  plastic 
cap  at  the  top.  A  back  pressure  of  0.5  bar  was  applied.  The  cell  pressures  are 
in  Table  2.3. 


2.6  TRIAXIAL  DEFORMATION 

For  the  drained  samples,  the  strain  rate  was  chosen  at  approximately  5  %  per 
day.  so  that  the  normally  consolidated  samples  would  fail  after  4  days  shear¬ 
ing,  approximately.  Despite  this  slow  strain  rate,  the  pore  pressures  some¬ 
times  rose  to  15  kPa  at  4  %  strain  approximately,  falling  thereafter. 

For  the  undrained  samples,  a  strain  rate  of  0.15  mm/min  was  chosen  to  limit 
the  effects  of  time  lag  in  the  pore  pressure  measuring  system,  equalization  of 
non-uniform  pore  pressures,  and  rate  dependent  effects. 


2.7  RECOVERY  OF  THE  SAMPLES 

The  measurements  of  the  microstructure  were  destructive,  so  it  was  necessary 
to  stop  the  triaxial  deformation  at  various  strains,  in  order  that  these 
measurements  could  be  made.  Control  samples  at  0  %  strain  were  obtained  by 
setting  up  samples  as  if  they  were  to  be  sheared,  omitting  all  shear  deforma¬ 
tion.  and  recovering  the  samples  immediately  after  the  triaxial  consolidation 
stage . 

After  stopping  the  triaxial  deformation  and  releasing  the  axial  load,  the  cell 
and  pore  pressures  were  reduced  slowly,  after  which  the  samples  were  removed 
from  the  triaxial  machines.  They  were  then  sliced  longitudinally  in  half  with 
a  cheesewire  along  the  plane  of  symmetry  if  one  had  formed.  One  half  of  the 
sample  was  used  to  determine  the  moisture  content,  and  the  other  half  was 
prepared  for  microstructural  examination  as  described  in  the  next  Chapter. 


oooooOooooo- 
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Fig.  2.1  Medium  pressure  filter  press. 
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Table  2.1  properties  01  Speswnite  kaolin 


Drained 

L’ndrained 

Pure 

S102 

47.0 

J6.2 

46 . 6 

by  wt 

A  i  2O3 

37.8 

38.  ~ 

39.5 

'' 

by  -At 

Fe:03 

0.6 

0.56 

0.0 

•J 

bv  wt 

T1O2 

0.09 

a 

by  wt 

Cao 

0. 20 

0 

by  wt 

MgO 

0.20 

■yr 

0 

by  wt 

k20 

)  1  .2 

1.01 

'0 

by  wt 

Na20 

) 

0.07 

by  wt 

H2O* 

13.1 

13.14 

14.0 

■y* 

by  wt 

PH 

5.0 

finer  than  10  pm 

99.5 

0 

by  wt 

finer  than  2  pm 

78-83 

■yf 

0 

by  wt 

mean  grain  size 

0."7 

pm 

specific  gravity 

2.6 

2.62 

liquid  limit 

69 

67.5 

mt 

plastic  limit 

34 

37.0 

or 

0 

plasticity  index 

35 

32.5 

•*r 

refractive  index 

1.56 

1.56 

*  or  loss  on  ignition. 


Table  2.2  Pressures  applied,  low  pressure  filter  presses,  bar 


_  Nominal 

Press  1 

Press  2 

Press  3 

hours 

®  0.125 

0.245 

0.245 

0.225 

24 

0.250 

0.370 

0.370 

0.350 

24 

1  0.500 

0.620 

0.620 

0.600 

24 

I  1.000 

1.120 

1.120 

1.100 

168 

^  Table  2.3  Cell  pressures  applied, 

,  undrained 

samples 

1  pressure,  bar 

1 

2 

4 

7 

®  duration,  day 

1 

1 

1 

2 

chapter  3 


ELECTRON  MICROSCOPY 


3.1  INTRODUCTION 

This  Chapter  reports  the  methods  of  electron  microscopy  and  of  image  analysis 
which  were  used  to  study  the  structures  of  the  samples  which  had  been  deformed 
in  the  triaxial  tests.  Virtually  all  of  the  electron  microscope  observations 
were  made  using  the  back-scattered  mode,  which  was  itself  relatively  new,  so 
some  improvements  in  preparation  techniques  had  to  be  made  for  it.  The 
automatic  methods  of  image  analysis  which  were  used  to  extract  quantitative 
information  from  the  electron  micrographs  were  almost  entirely  novel:  and  the 
development  and  testing  of  these  new  methods  was  a  major  task.  These  methods 
are  discussed  in  turn  below:  Fig.  3.1  shows  a  typical  scanning  electron  micro¬ 
graph;  and  Figs.  j. 2  and  3.3  summarize  the  principal  analyses  made  in  Glasgow 
and  UEA,  respectively. 


3.2  SCANNING  ELECTRON  MICROSCOPY 
3.2.1  Introduction 

The  first  steps  are  to  prepare  the  samples  for  scanning  electron  microscopy, 
to  examine  them  in  the  microscope,  and  to  digitize  the  electron  micrographs. 
This  Section  describes  these  steps. 


3.2.2  Impreenat ion 

At  tht  end  of  each  triaxial  test,  the  samples  were  first  sliced  into  two 
halves  ior.gitudinal ly,  as  described  in  the  last  Chapter.  Then,  the  samples 
were  impregnated  following  the  methods  described  by  Smart  and  Tovey  (1982). 

The  usual  procedure  for  the  drained  tests  was  to  cut  sub-samples  sized  20  mm 
square  by  10  ram  thick.  These  were  placed  face  down  in  a  foil  dish  and  immersed 
directly  in  acetone.  The  acetone  was  changed  several  times  at  24  hour  inter¬ 
vals.  When  all  water  had  been  removed  by  the  acetone,  a  mixture  of  Araldite 
AY18 ,  Hardener  HZ18,  and  1-2  %  di-n-butyl  phthalate,  diluted  with  an  equal 
volume  of  acetone,  was  introduced  into  the  sample.  The  acetone  was  allowed  to 
evaporate  slowly  from  the  mixture,  which  was  topped  up  as  necessary  to  keep 
the  sample  submerged;  then  the  resin  was  left  for  2-3  days  to  harden.  Finally, 
it  was  cured  at  40  degrees  C. 

For  the  undrained  tests,  complete  half  samples  were  usually  impregnated.  These 
were  placed  face  down  in  a  small  dish  and  packed  round  with  sand;  the  small 
dish  was  then  placed  in  a  larger  dish  litted  with  a  drain.  Both  dishes  were 
carefully  flooded  with  25  %  acetone  in  water;  then,  at  intervals  of  one  or  two 
days,  the  outer  dish  was  drained  and  refilled  in  turn  with  50%,  75  %  and  100 


%  acetone.  At  least  three  changes  of  100  %  acetone  were  used.  When  ail  the 
water  had  been  removed  by  the  acetone,  a  mixture  of  100  parts  Vestopal  W.  1 
part  initiator,  0.5  parts  activator,  (by  weight),  and  200-300  parts  acetone 
(by  volume),  was  added  to  the  small  dish.  The  acetone  was  allowed  to  evaporate 
slowly  from  the  mixture,  which  was  topped  up  as  necessary  to  keep  the  sample 
submerged;  then  the  resin  was  hardened  at  room  temperature.  Finally,  it  was 
cured  at  60 *C  for  48  hours. 

Especially  during  the  early  stages  of  the  study,  ethanol  and  methanol  were 
used  to  remove  the  pore  water,  and  then  acetone  was  used  to  replace  the 
methanol . 


3.2.3  Grinding 

When  the  impregnated  samples  had  hardened,  2  mm  was  ground  off  the  face  to 
remove  the  disturbed  zone.  Then  sub-samples  were  cut,  ground  and  polished  to 
prepare  a  block  for  observation  in  the  scanning  electron  microscope.  The  final 
surface  should  be  flat  and  free  not  only  from  scratches  and  pluck  marks  but 
also  from  smearing.  Conventional  methods  of  grinding  thin  sections  for  optical 
microscopy  do  often  leave  a  thin  smeared  zone  at  the  surface;  and  the  require¬ 
ments  for  scanning  electron  microscopy  are  more  demanding.  The  success  of  this 
project  owes  much  to  the  skill  and  patience  of  Mr . I .Marshal  1 .  who  was 
responsible  for  improving  the  final  preparation  technique  and  preparing  the 
blocks?  The  samples  were  ground  successively  on  emery  cloth  Nos.  240  and  600 
lubricated  by  water  and  6  pm  and  1pm  diamond  lubricated  with  Hyprez  fluid 
(supplied  by  Engis  Ltd.). 


3.2.4  Coating 

The  polished  blocks  were  mounted  on  specimen  stubs  using  silver  dag  and  then 
coated  with  carbon. 


3.2.5  Observation 

For  the  drained  samples,  almost  all  of  the  scanning  electron  microscopy  was 
done  using  the  Hitachi  L-4  0  SEM  in  the  School  of  Environmental  Sciences  in 
UEA.  Digitized  images  were  acquired  by  a  Link  Analytical  interface  (which  is 
also  capable  of  X-ray  micro-anal}’sis 1 .  converted  to  Semper  format  by  a  PC.  and 
fed  into  the  image  archiving  system  described  below.  The  whole  route  from  SEM 
through  to  the  archive  is  hard-wired.  For  the  undrained  samples,  electron 
micrographs  were  digitized  and  recorded  by  a  computer  which  is  an  integral 
part  of  the  Cambridge  360  SEM  in  Glasgow;  they  were  then  passed  to  a  second 
computer  which  is  close-coupled  to  the  first,  where  they  were  stored 
temporarily  on  hard  disk.  At  the  end  of  a  session,  the  micrographs  were  down¬ 
loaded  to  a  half-inch  magnetic  tape,  which  can  held  about  250  micrographs  and 
forms  the  primary  archive.  Duplicate  copies  of  about  5%  of  the  images  were 
stored  on  floppy  discs  for  the  Nimbus  image  analyzer. 


The  main  series  of  micrographs  were  obtained  using  the  back-scattered  mode; 
some  additional  micrographs  were  obtained  using  the  secondary  electron  mode 
for  comparative  purposes.  Untilted  samples  were  used.  This  technique  presents 


*  Mr.  A.  A.  Burnett  prepared  the  blocks  in  Glasgow. 
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a  cross-sectional  view  of  the  sample,  vnich  simplifies  the  image  analysis, 
'ontrast  in  the  bacx-scat tered  mode  arises  primarily  from  atomic  number  con¬ 
trast.  so  that  particles  in  general  appear  bright  against  a  dark  background  of 
resin-filled  pores.  Fig. -3.1.  The  incident  electrons  penetrate  some  distance 
mto  the  sample,  and  there  is  some  spreadine  of  electrons  within  the  sample 
before  they  escape  back  to  the  collector.  Thus,  some  particles  which  are 
actually  below  the  true  surface  may  lie  seen  more  or  less  faintly:  and  also 
there  is  some  softening  of  the  contrast. 


All  the  digitized  micrographs  (images)  taken  for  this  study  have  512  rows  of 
512  pixels  each  and  are  digitized  to  256  gray  levels.  Usually  the  images  were 
taken  from  a  vertical  face  aligned  as  far  as  possible  with  a  plane  of  symmetry 
in  the  sample.  In  UEA.  34  images  were  usually  taken  from  each  sample  as  fol¬ 
lows  : 


number 
of  images 
24 
6 
4 


nominal 
magnificat  ion 
x2000 
xlOOO 
x500 


image  width 
mm 

0.056 

0.113 

0.225 


In  Glasgow,  26  images  were  usually  taken  from  each  sample  as  follows: 


number 
of  images 
25 
l 


nominal 
magnificat  ion 
X4000 
xSOO 


image  width 
mm 


0.025 


0.120 


The  main  analyses  were  made  using  the  x2000  and  x4000  images,  with  the  others 
held  in  reserve.  As  far  as  possible,  the  images  in  UEA  were  spaced  at  1  mm 
centers  in  two  traverses  of  12  images,  one  being  horizontal  and  the  other 
vertical,  the  images  in  Glasgow  were  spaced  at  2  mm  centers  in  a  5  x  5  array. 
In  UEA,  with  the  exception  of  some  of  the  earlier  samples,  the  contrast  was 
standardized  at  the  start  of  each  batch  of  images  by  observing  a  specially 
prepared  block,  which  contained  large  areas  of  clear  Araldite  and  large  quartz 
grains.  The  gray  level  on  the  Araldite  was  adjusted  to  0.  and  that  on  the 
quartz  to  221.  Images  of  the  areas  used  for  this  adjustment  were  also 
digitized. 


3.3  IMAGE  ANALYSIS 
3.3.1  Introduct ion 

At  the  start  of  this  project,  both  universities  installed  PC-based  versions  of 
the  Semper  image  analysis  system;  for  convenience,  the  PC  in  Glasgow  is  called 
Nimbus.  These  effectively  replaced  the  main-frame  and  Quantimet  systems  which 
had  been  in  use  previously.  By  now,  the  original  installations  have  been  pro¬ 
gressively  improved.  In  Glasgow,  additional  support  from  SERC  led  to  the  use 
of  transputers.  In  UEA,  net-worked  PCs  and  a  Sun  workstation  have  been  added, 
all  running  Semper. 
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3.3.2  Image  archives 

Each  digitized  image  required  0.25  MB  storage  space,  in  round  terms,  each  oi 
the  three  major  experiments  used  40  samples  and  25  images  from  each.  i.e.  1000 
SEM  images  per  experiment,  excluding  all  subsidiary  and  supporting  images  and 
experiments.  Each  of  these  experiments  also  used  about  1000  optical  images.  In 
addition,  there  has  been  a  necessity  to  store  processed  images  ready  for  sub¬ 
sequent  stages  of  analysis:  and  the  data  bank  of  raw  results  is  also  large. 

The  image  analysis  facility  in  UEA  by  now  consists  of  an  network  of  S02S6. 
80386.  and  80486  PCs.  all  connected  by  D-link.  so  that  each  computer  can  read 
the  hard  disk  of  any  other.  Fig.  3.4.  The  network  also  contains  a  Sun  worksta¬ 
tion.  a  Gigabyte  hard  disk,  optical  disk  reader,  and  tape  drive.  The  SEW 
itself  is  connected  to  the  net  as  indicated  above:  and  there  is  another  con¬ 
nection  to  the  Cn i vers i tv’s  main  network  and  central  Vax  computer.  Thence, 
there  is  a  JANET  connection  to  Glasgow.  These  facilities  have  been  expanded 
gradually,  and  there  have  been  some  changes  in  practice  as  new  facilities  were 
added.  In  general,  images  are  kept  on  the  Gigabyte  hard  disk  and  copied  from 
there  to  the  computers  when  required.  Duplicate  copies  are  kept  in  the  archiv¬ 
ing  system  connected  to  the  Vax:  and  third  copies  of  a  large  proportion  of  the 
collection  have  been  stored  on  optical  disks  or  magnetic  tapes.  The  archiving 
system  has  not  been  trouble  free,  and  from  our  experience,  we  are  inclined  to 
recommend  that  important  data  should  be  held  in  triplicate.  All  the  images 
originating  from  UEA  are  stored  in  Semper  format. 

Early  in  the  project,  some  images  were  transferred  from  UEA  to  Glasgow  by 
electronic  transfer  from  UEA’s  Vax  to  Glasgow’s  ICL  mainframe.  However,  for 
small  quantities  of  representative  or  special  images,  it  was  more  convenient 
to  send  floppy  disks  by  post;  and  for  large  quantities,  images  were  copied 
from  UEA’s  Vax,  onto  Vax-type  tapes,  sent  to  Glasgow,  read  into  a  Vax  there, 
and  transferred  thence  to  the  Meiko  computer.  Images  sent  from  UEA  are  stored 
in  home/psmart/uea. 

In  Glasgow,  when  each  tape  from  the  SEM  is  full,  the  Computing  Service  copies 
it  via  a  Sun  workstation  and  Ethernet  to  the  CPU  of  the  Meiko  (which  is  an 
embedded  Sun  board),  and  thence  to  the  Meiko’s  main  disk,  where  they  are 
stored  in  directories  called  home/psmart/bai .  psmart  is  used  for  archiving 
only,  and  all  workers  in  the  group  have  read-only  access  to  it.  A  separate 
sub-directory,  /tapeN.  where  N  is  a  single-digit  number,  holds  all  the  images 
from  one  half-inch  tape.  Due  to  peculiarities  in  the  file-handling  systems, 
the  Glaswegian  micrographs  were  reordered  before  being  down-loaded  to  the 
half-inch  tapes,  and  the  original  file  name  was  moved  into  the  file;  there¬ 
fore.  files  are  given  new  names,  /fileDDD,  where  DDD  is  a  three  digit  number 
starting  at  001.  The  format  then  becomes:  header  (file  name)  512  bytes;  header 
(original)  512  bytes;  512  lines  x  512  bytes  each.  All  the  files  from  one  tape 
are  then  inspected:  the  images  are  presented  on  one  screen  at  about  4  second 
intervals,  whilst  the  new  and  original  file  names  are  shown  on  another  screen 
and  listed  to  file.  The  list  enables  the  images  to  be  reordered  correctly  by 
logging  onto  psmart  and  using  a  home-made  Unix  command  to  swap  pairs  of  file 
numbers;  and  the  visual  presentation  enables  bad  images  to  be  spotted  and 
moved  to  the  ends  of  the  batches  containing  them.  After  the  final  reordering, 
the  scanning  program  is  used  again  to  list  the  names  and  numbers  of  the  images 
to  a  catalogue;  comments  are  edited  in;  and  the  catalogue  is  saved.  During 
analysis  from  a  Sun  workstation,  two  windows  can  be  opened,  both  can  log  onto 
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different  sets  of  transputers,  me  of  wnich  cun  ao  the  work.  Ahiist  the  other 
reaas  the  catalogue.  Fie. -.''.5. 


3.3. 3  Image  analysis  hardware  and  software 

In  UEA,  the  image  analysis  hardware  consists  of  the  net-worked  computers 
described  above.  Most  of  the  PCs  are  equipped  with  frame  stores  and  secondary 
monitors  to  display  the  images:  the  Semper  software  image  processing  routines 
can  be  run  on  any  of  these  PCs.  The  Sun  workstation  has  its  own  Semper  soft¬ 
ware  and  uses  a  window  on  its  main  monitor  for  displaying  images.  There  is  a 
considerable  difference  in  processing  speed  between  the  fastest  and  the 
slowest  of  these  computers:  and  this  is  taken  into  account  when  allocating 
vork  to  them.  The  Sun  is  capable  of  running  several  processes  in  parallel:  and 
it  has  generally  proved  to  be  faster  (in  terms  which  include  operator  time)  to 
place  several  medium  sized  standard  jobs  on  it  rather  than  one  very  large  one. 

All  image  analysis  in  UEA  is  done  using  the  Semper  system.  Apart  from  interac¬ 
tive  use  for  occasional  tasks,  programs  are  written  in  the  Semper  language  and 
left  to  run  unattended,  perhaps  for  24  hours  or  overnight.  There  are  facili¬ 
ties  in  the  Semper  system  both  to  call  Fortran  sub-routines  which  are  external 
to  the  Semper  package  proper  nr  to  import  these  sub-routines  into  the  Semper 
package  as  extra  ‘verbs’  (commands).  The  work  has  involved  a  combination  of 
all  these  methods  of  operation. 

A  different  course  of  development  was  followed  at  Glasgow.  During  the  first 
year  of  the  project,  an  80386  PC  with  Semper  frame-store  was  installed.  This 
was  used  for  display;  but,  as  a  matter  of  convenience,  external  Fortran 
programs  were  used  for  the  novel  algorithms  which  were  being  developed.  These 
new  algorithms  were  running  very  slowly;  and  there  was  a  memory  restriction 
even  though  the  code  had  been  written  to  read-and-process  line-by-line.  At 
this  stage,  SERC  loaned  another  PC  with  a  T4  transputer,  and  the  improved 
speed  and  simplified  management  of  memory  provided  an  excellent  development 
environment.  This  was  upgraded,  first  to  4  T4  transputers,  then  to  6  T8 
transputers  with  a  second  monitor,  and  finally  to  the  Meiko  computer  (which  is 
a  University  central  facility).  Full  advantage  has  not  yet  been  taken  of  the 
Meiko:  because  the  main  programs  written  for  the  6-transputer  system  were  fast 
enough,  and  so  priority  has  been  given  to  algorithm  development.  In  practice, 
although  the  IPLIB  image  analysis  package  is  available,  all  code  used  for  the 
Transputers  has  been  specially  written,  because  so  much  of  it  was  novel. 

The  general  arrangement  of  (our  part  of)  the  Meiko  is  shown  in  Fig. -3.5.  The 
’top’  or  ‘master’  transputer  communicates  with  the  terminal,  reads  images  from 
disk,  passes  these  and  commands  to  the  other  transputers,  receives  results, 
and  combines  these  results  as  necessary.  The  four  central  transputers, 
’workers’,  each  receive  one-quarter  of  the  image  and  process  that;  inter¬ 
mediate  results  are  swapped  between  these  processors  as  necessary.  The  ’dis¬ 
play’  transputer  presents  images  on  the  secondary  monitor;  it  may  also  do  some 
processing  of  the  recombined  computed  images.  Each  transputer  has  4  MB  RAM, 
which  is  barely  enough.  At  one  time,  separate  versions  of  the  programs  were 
maintained  to  analyze  either  one  or  a  series  of  micrographs;  but  most  of  the 
single-image  programs  are  being  phased  out.  Future  programs  will  probably  con¬ 
tain  pre-processing  transputers  before  the  workers,  and  post-processing 
transputers  after  them;  although  this  will  increase  the  starting  up  time. 
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iodine  tnese  extra  transducers  into  tne  uipe  w  i  i  i  not  affect  the  rate  of 
process ins.  Fhe  Meiko  is  now  usually  run  irom  a  run  works tat ion.  so  it  is  pos¬ 
sible  to  ioe  onto  two  sets  of  transputers  at  once,  which  helps  when  deveiopine 
parallel  programs,  reading  the  file  catalogue,  etc. 

After  returning  SERC’s  PC  and  transputers,  .‘limbus  was  fitted  with  a  single  I'S 
transputer  and  a  third  monitor,  see  Fig.  3.5.  At  present,  there  is  a  choice  of 
using  Semper  or  tiie  transputer:  in  theory,  by  using  windows,  both  systems 
could  be  used  at  once.  New  programs  are  usually  written  on  the  transputer,  for 
vhich  there  are  Occam.  Fortran,  and  C  compilers:  but  most  work  uses  both 
systems.  All  of  the  optical  microscopy  reported  in  the  next  chapter  used  this 
PC.  Images  were  digitized  by  Semper  using  a  thermionic  camera,  which  was 
preferred  to  avoid  the  problems  of  re-digitizing  which  arise  unless  special 
precaut ions  are  taken  with  CCD  cameras:  and  the  calculations  were  made  by  the 
transputer.  These  images  were  stored  on  floppy  disks  after  the  calculations 
had  been  completed. 


3.4  IMAGE  RECTIFICATION 
3.4.1  introduct ion 

Before  analyzing  the  digitized  electron  micrographs,  a  number  of  steps  were 
taken  to  rectify  defects  in  them.  These  steps  are  discussed  in  turn  below. 


3.4.2  Impregnation  and  grinding  defects 

Occasionally,  during  scanning  electron  microscopy,  the  location  of  the 
recorded  area  was  changed  slightly  to  avoid  patches  where  either  the  impregna¬ 
tion  or  the  grinding  was  defective. 


3-4.3  Charging 

Occasionally,  during  scanning  electron  microscopy,  some  areas  of  the  samples 
charged  up  and  became  brighter  than  they  should  have  been.  These  areas  were 
few  and  small;  moreover,  both  of  the  principal  methods  of  analysis  used  here, 
i.e.  intensity  gradient  methods  and  the  method  used  to  obtain  black-and-white 
images,  minimize  the  effect  of  charging.  Therefore,  no  special  measures  were 
employed. 


3.4.4  Mineralogical  variations 

A  few  of  the  clay  particles  appeared  brighter  than  average:  this  was 
attributed  to  differences  in  composition  such  as  an  increase  in  iron  content. 
The  effect  was  small  and  rare  and  no  special  measures  were  employed. 
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3.4.5  ■  jver-run 

occasionally,  some  of  the  particles  looked  as  if  they  nad  very  small  holes 
vvithin  them.  This  was  attributed  to  malfunction  of  the  analouue-to-digitai 
converter:  i.e.  pixels  which  should  have  been  i maxsrev  +  ru  were  recorded  as 
n.  In  other  words,  instead  of  being  digitized  to  8  bits,  the  intensities  had 
been  digitized  to  d  bits,  and  the  most  significant  bit  had  been  discarded.  A 
short  Semper  program  was  written  to  patch  these  holes  empirically:  and  a  sec¬ 
ond  program  was  used  to  check  whether  this  had  been  successful. 


3.4.6  Electrical  interference 

During  the  course  of  the  investigation,  it  became  apparent  that  the  electrical 
circuits  within  the  scanning  electron  microscope  and  digitizer  were  being 
affected  by  electrical  interference,  either  from  within  the  instruments,  or 
from  nearby  apparatus,  or  from  the  power  supplies.  This  interference  was  dis¬ 
covered  when  Fourier  transforms  of  images  of  plain  resin  blocks  were 
inspected:  the  power  spectrum  was  found  to  have  large  peaks,  which  were  un¬ 
related  to  the  structureless  blocks,  and  which  usually  occurred  in  the  same 
places.  These  peaks  were  also  found  in  the  power  spectra  of  the  micrographs  of 
the  samples  proper.  Similar  effects  have  since  been  reported  to  us  by  other 
image  analysts. 

To  mitigate  this  artifact,  the  power  spectrum  was  searched,  and  the  peaks  were 
located  and  masked  out.  The  rectified  image  was  then  recovered  by  taking  the 
inverse  Fourier  transform.  In  practice,  the  power  spectrum  has  a  natural  peak 
at  its  center;  so  the  masking  was  confined  to  the  outer  part  of  the  power 
spectrum. 

This  operation  was  combined  with  a  deconvolution  of  the  point  spread  function, 
see  wiener  filtering  below. 


3.4.7  Signal-to-noise  ratio 

Close  inspection  of  the  electron  micrographs  showed  that  they  were  granulated 
very  finely  due  to  high  frequency  noise.  This  noise  was  measured  for  each 
batch  of  micrographs  by  finding  the  variance,  ,  of  the  gray  levels  in  clear 
areas  of  the  control  micrographs  which  had  been  taken  on  the  plain  resin 
blocks. 

To  obtain  the  signal-to-noise  ratio  of  an  ordinary  micrograph,  the  variance, 
vt ,  of  its  gray  levels  was  found.  Thence,  the  signal  to  noise  ratio,  sT ,  was 
calculated  from: 


Sr  =  sqrt(  (vt  -  vh)  /  vh  )  (3.1). 


3.4.8  Point  spread  function 

When  scanning  electron  micrographs  are  being  formed,  an  electron  beam  is 
directed  at  the  sample.  Some  of  the  incident  electrons  are  reflected  as  back- 
scattered  electrons;  some  liberate  secondary  electrons  of  low  energy;  and  some 
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ire  absorbed,  in  che  back-:>cat tereu  mode  of  operation,  both  the  back-scattered 
and  the  secondary  electrons  are  collected  and  contribute  to  the  intensity. 
Although  the  incident  beam  is  of  very  small  diameter,  the  electrons  in  it 
enter  the  sample  and  are  spread:  so  the  received  signal  contains  contr lbut ions 
from  regions  around  the  theoretical  point  of  contact  (Smart  and  Tovey,  19S2. 
Chapter  M.  This  effect  leads  to  a  blurring  of  the  image  and  some  loss  of  con¬ 
trast.  The  mathematical  description  of  the  blurring  is  known  as  the  optical 
transfer  function  or  as  the  point  spread  function.  For  an  electron  microscope, 
it  is  expected  that  the  shape  is  Gaussian  (Lewis  and  Sakrison,  19^5:  Guan  and 
Ward.  1989). 

The  point  spread  function  was  measured  from  electron  micrographs  of  small  par¬ 
ticles  of  dirt  which  were  seen  on  the  plain  resin  blocks.  The  particles 
selected  were  small  enough  to  be  effectively  point  sources.  Their  intensity 
profiles  were  measured,  averaged  rotational ly  about  the  peak  positions,  and 
then  averaged  collectively.  A  Gaussian  profile  was  fitted  to  this  final 
average  profile.  This  gave  an  overall  spread  diameter  of  11.4  pixels  at  x2000 
magnification,  corresponding  to  a  Gaussian  rms  radius  of  11.4/8  pixels  or 
0. 156  microns. 


3.4.9  Wiener  filtering 

After  masking  the  electrical  interference  from  the  Fourier  transform  of  an 
image,  Wiener  filtering  was  used  to  sharpen  the  image  by  removing  the  blur: 

F(u.i')  =  G(  u,  v)  H*  (  u,  v)  /  {mod(H(  u.  v)2  )  +  l/sr2}  (3.2). 

where  G(u.v)  is  the  masked  power  spectrum  of  the  image; 

H( u. v)  is  the  power  spectrum  of  the  point  spread  function; 

H *{u.v)  is  the  complex  conjugate  of  the  point  spread  function. 

The  rectified  image  was  recovered  from  F(u.v)  by  taking  the  inverse  transform. 
An  example  is  shown  in  Fig.3.7.  Fast  Fourier  transforms  were  used  throughout 
these  calculations. 

(Further  details  are  given  by:  Stroke  et  ai..  1971.  who  used  optical  process¬ 
ing:  Lewis  and  Sakrison,  1975;  Bates  and  McDonnel,  1986;  Gonzalez  and  Wintz, 
1987;  Guan  and  Ward,  1989. ) 


3.4.10  Missing  scan  lines 

Some  of  the  electron  micrographs  taken  in  Glasgow  were  found  to  contain  miss¬ 
ing  scan  lines.  During  the  inspection  mentioned  above,  for  each  image,  the  sum 
of  the  gray  levels  is  calculated  for  each  line  of  the  file  including  the  two 
headers.  Confidence  limits  for  these  sums  are  set  at  three  standard  deviations 
from  the  mean.  In  our  case,  the  headers  should  lie  outside  the  limits;  so  the 
program  shows  (a)  how  many  header  lines  lie  within  the  limits;  and  (b)  how 
many  ordinary  lines  lie  outside  them.  Suspect  images  are  re-examined  one-by- 
one  and  removed  if  necessary.  At  the  time  of  writing,  some  of  these  images 
have  also  been  found  to  have  'jumbled’  scan  lines,  which  are  found  during  the 
visual  examination;  and  we  are  considering  using  similar  confidence  limits  for 
’mean  chord  lengths’  to  search  for  these. 
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.5  INTENSITY  GRADIENT  ANALYSIS 

The  clay  used  in  these  studies  consists  a i most  entire!:'  of  flat  plate- .  The 
electron  micrographs  show  a  planar  cross-section  through  the  material,  in 
which  most  of  the  plates  appear  as  thick  lines,  .-.part  from  a  few  areas  .^here 
the  piates  are  jumbled  together  randomly,  there  is  generally  strong  iocai 
preferred  orientation:  many  of  the  micrographs  also  show  an  overall  preferred 
orientation.  Intensity  gradient  analysis  was  used  to  measure  both  the  strength 
and  direction  of  this  overall  preferred  orientation. 

In  an  electron  micrograph  showing  bright  particles  against  a  dark  background, 
the  intensity  gradient,  grad  /.  is  a  vector,  whose  value  is  high  at  the  edges 
of  the  particles,  and  whose  direction  s  perpendicular  to  those  edges.  The 
analysis  starts  by  calculating  the  intensity  gradient: 


u  = 

grad  / 

<3. 

.3) 

V  - 

mod  grad  I 

(3. 

.4) 

.4  = 

arg  grad  1 

(3, 

.5) 

As  a  general  rule,  the  20,14-formula  is  used,  because  this  formula  appears  to 
suit  the  work  in  hand.  The  formula  uses  the  20  nearest  pixels  and  was  derived 
from  a  least  squares  fit  for  the  first  14  differential  coefficients  in  a  two 
dimensional  Taylor  series  (Smart  and  Tovey,  1988a).  To  avoid  cases  close  to  A 
=  atan  (0/0).  pixels  are  labelled  ’undecided’  if  U  is  smaller  than  an 
arbitrary  limit.  The  default  value  for  this  limit  is  3  pixels,  where,  for  con¬ 
venience,  the  term  ’pixel’  is  used  as  a  measure  of  length:  but  higher  values 
may  sometimes  be  needed. 

The  results  are  often  summarized  in  an  unweighted  rosette  frequency  histogram 
plot  of  orientation  for  each  image.  Fig. -3. 8. a.  Earlier  studies  using  this 
approach  often  gave  spiky  rosettes:  but  a  combination  of  a  better  digitiza¬ 
tion.  a  more  sophisticated  formula,  better  choice  of  the  limit  for  the 
undecided  case,  and  a  more  appropriate  histogram  bin  size,  has  resulted  in 
reasonably  smooth  rosettes.  The  best  fitting  ellipse  was  found  for  the  rosette 
for  each  image,  and  an  Index  of  Anisotropy,  /A ,  found  from: 

/a  =  (major  axis  -  minor  axis)/  minor  axis  (3.6). 


3.6  MAPPING  BY  TOP-CONTCURING 

The  clay  structure  consisted  mainly  of  groups  of  particles  showing  strong 
local  preferred  orientation  which  are  called  domains  here.  There  were  also 
some  random  clusters  of  plates.  Two  new  algorithms  were  developed  to  map  these 
features  ready  for  measurement , 

The  first  of  these  algorithms  to  be  developed  maps  domains  and  random  clusters 
by  first  encoding  A  into  a  limited  number  of  directions,  and  then  smoothing 
the  encoded  image  (Smart  et  al.  1988).  Typically,  4,  8,  or  12  directions  are 
used.  The  first  successful  implementation  used  a  T414  transputer  in  a  PC:  but 
the  main  development  used  six  T800  transputers  in  a  PC.  This  algorithm  has  now 
been  moved  to  the  Meiko  without  change  apart  from  some  simplification  and  the 
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use  or'  a  different  erupnics  board.  A  Semper  version  nas  been  -Aril  ten:  and  this 
ls  now  available  via  Synoptics  Ltd. 

The  encoded  image  is  smoothed  by  selecting  the  most  popular  direction  from  the 
histoeram  of  encoded  angles  within  a  large  circular  uniform  filter  (see  Appen¬ 
dix  3.1).  The  filter  is  topically  of  20  pixels  ^adius.  Reflection  is  used  to 
avoid  leaving  unprocessed  borders.  If  more  than  an  arbitrary  number  of  pixels 
within  the  filter  have  undecided  angles,  then  the  pixel  in  the  map  is  labelled 
undecided;  by  default,  the  code  uses  the  radius  of  the  filter  as  this  limit, 
and.  m  the  current  work,  it  is  rarely  if  ever  exceeded.  The  pi. .el  in  the  map 
is  labelled  'random'  if  the  difference  between  the  maximum  and  the  mean  of  the 
histogram  is  below  a  limit.  EXC  wTu^-h  may  be  set  arbitrarily,  but  which  .s 
by  default  set  by: 

EXC  =  200.* SQRT ( -LOG ( PROB ) /X I F ) * NOD* S I N ( P I /NOD )  /  PI  (3.7). 

where  MOD  is  the  number  of  directions.  NIE  is  the  number  of  pixels  in  the  fil¬ 
ter,  and  PROB  is  the  probability. 

The  difference  between  the  maximum  and  the  mean  is  thought  to  be  more  robust 
than  the  difference  between  the  maximum  and  the  minimum  in  cases  in  which  only 
one  direction  differs  from  the  others.  A  bias  may  arise  if  the  histogram  has 
two  or  more  equal  maxima,  and  *he  present  code  contains  one  of  several  algo¬ 
rithms  which  were  considered  to  minimize  this  bias. 

As  far  as  we  know,  no  similar  jnalysis  had  ever  been  reported;  so  some  time 
was  spent  in  testing  the  values  of  the  various  limits,  in  deciding  the  radius 
of  the  filter  and  how  many  directions  to  use,  and  in  deciding  the  magnifica¬ 
tion  and  the  other  instrumental  settings  on  the  microscope. 

The  most  difficult  parameter  to  fix  was  the  radius  of  the  filter,  because 
there  was  originally  thought  to  be  no  predetermined  optimum  size  for  it.  The 
only  rule  found  was:  v.o  smooth  away  features  of  a  given  size,  use  a  filter  of 
a  slightly  larger  size.  In  the  main  series  of  micrographs  to  be  analyzed, 
there  were  features  of  all  sizes,  so  the  final  choice  of  radius  for  the  filter 
was  made  on  the  basis  of  visual  assessment  by  a  p>nel  of  experienced  micro- 
scopists.  The  nominal  value  chosen  was  19  pixels:  but  the  code  increases  the 
nominal  radius  by  0.4  pixels  to  improve  the  circularity  of  filters  less  than  8 
pixels  radius.  Subsequently,  a  modified  chi-square  test  was  used  to  test  this 
choice.  The  image  was  mapped  at  various  radii;  the  areas  were  averaged  over 
all  radii;  and  the  sum  of  the  squares  of  the  errors  of  the  original  areas  from 
the  corresponding  averages  was  treated  as  an  error  function.  The  result  for  a 
typical  image  is  in  Fig. -3. 9. a.  The  choice  of  19  p'xels  was  also  supported  by 
some  consistency  ratio  results  and  by  a  preliminary  study  of  the  semi- 
variogram.  both  of  which  are  discussed  below. 

This  algorithm,  which  was  called  topcontour ing  by  its  author,  was  a  great  step 
forward  in  the  field  of  soil  microstructural  analysis.  It  is  faster  than  the 
alternative  method  described  in  the  next  Section. 
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CONSISTENCY  RATIO  MAPPING 

' . " . 1  ■  onsistencv  ratio  mapping  per  ce 

Consistency  ratio  mapping  is  the  second  of  the  new  algorithms  to  he  developed 
for  mapping  clay  structures.  It  is  more  accurate  than  topcontour ing.  .-slightly 
slower,  and  requires  more  memory:  it  maps  domains  and  random  clusters  by  first 
smoothing  the  vector  field.  U.  then  testing  for  randomness  and  segmenting  the 
oriented  areas  into  a  limited  number  of  equally  spaced  directions  (Smart. 

14QQ.  Smart  &-  Leng.  I*J91).  .An  example  is  shown  in  Fie. 3.1.  The  equations  are: 


.Y 

- 

U  cos  2.1 

(3.8) 

r  r 
t 

= 

r* 

U  sin  2.4 

(3.9) 

s 

= 

- 

U 

(3.10) 

K 

= 

V 

i  A2  +  r2  ) 

(3.11) 

c 

- 

R/S 

(3.12) 

T 

= 

1 

> 

a  tan  (  17  A) 

(3.13) 

C  is  called  the  consistency  ratio  (see  Smart  and  Tovey.  1982.  Section  12.1.4 
for  a  review).  If  5  is  small,  T  must  be  labelled  undecided  and  mapped  accor¬ 
ding  1  y .  I  f  : 


C  <  vr{  ( 1/AO  In  (1/P)  }  (3.14), 

where  .V  is  the  number  of  (decided)  pixels  in  the  filter,  and  P  is  the 
probability  of  obtaining  a  greater  consistency  ratio  by  chance,  then  the  pixel 
is  mapped  as  random. 

For  analyzing  clay  structure,  the  angle  is  multiplied  by  2  at  the  start  of  the 
calculation  and  divided  by  2  at  the  end.  because  the  small  linear  features 
have  two-fold  symmetry.  (Smart  and  Leng,  1991,  have  extended  the  concept  to 
chromatic  mapping,  which  lacks  this  symmetry.) 

In  addition  to  making  the  calculations  within  the  smoothing  filter,  they  were 
also  made  for  each  image  and  for  each  sample  to  provide  an  alternative  measure 
of  preferred  orientation. 

Further,  calculating  consistency  ratio  as  a  function  of  size  of  filter  showed 
that  it  fell  rapidly  at  first  as  the  size  of  the  filter  increased,  then  it 
appeared  to  level  off  at  a  level  representative  of  the  image  as  a  whole.  In 
Fig. 3. 9. b.  which  is  for  an  image  of  nominal  magnification  of  x500,  the  break¬ 
point  seems  to  occur  at  about  5  pixels  radius,  corresponding  to  20  pixels 
radius  for  a  magnification  of  x2000.  This  result  gave  some  further  support  for 
the  choice  of  size  of  filter  for  the  main  analysis  (see  also  Section  3.10.6 
below). 

Consistency  ratio  mapping  was  implemented  on  the  Meiko  by  rewriting  the  rele¬ 
vant  parts  of  the  kernel  of  the  topcontouring  code.  The  memory  requirement  is 
4  MB  per  transputer  (see  the  method  of  calculation  given  by  Smart  and  Leng, 
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I ->90a ) .  The  Semper  version  is  .1  direct  development  i  *. *:e-  .  cp-  e  r.teer  me 
procram. 

Although  consistency  ratio  mappine  .  •  ?  »e/.ver  than  top- contouring.  .  i  .s 
hemet  i till  iy  preferable  and  easier  to  set  up.  ince  less  parameters  -save  to 
•e  chosen.  The  main  analyses  reported  below  usen  consistency  rat  10  mapping 
'ith  a  r.inius  of  !’>  pixels.  The  results  presented  elow  are  i  or  domains  mapped 
*nto  4,  a  or  13  directions:  there  is  also  <i  random  c lass.  Although  some  points 
in  the  preceding  intensity  gradient  analysis  must  he  labelled  undecided,  the 
number  of  pixels  which  remained  undecided  after  either  top-contouring  or  con¬ 
sistency  ratio  mapping  was  either  zero  or  negligible.  Mapping  into  four  direc¬ 
tions  allows  the  original  images  to  be  colored:  e.g.: 


horizontal  domains 
down- left  domains 
vertical  domains 
down- right  domains 
random  clusters 


green 

blue 

red 

ye  1  low 

purple  or  white 


In  UEA,  the  consistency  ratio  maps  are  compressed  by  chord  length  coding  and 
stored:  the  practice  in  Glasgow  is  to  re-calculate  if  necessary. 


3 .  .  2  Post-processing 

After  the  electron  micrographs  had  been  mapped,  the  mapped  features  were 
measured  using  the  standard  routines  provided  by  the  Semper  image  analysis 
package.  To  assess  the  capability  of  this  package,  and  to  compare  alterna¬ 
tives,  25  separate  measurements  were  made  on  each  domain  of  each  image  for 
both  drained  series  (see  Appendix  3.2).  Many  of  these  measurements  were  effec¬ 
tively  alternatives;  for  example,  Feret  diameters,  equivalent  circular 
diameter,  and  mean  square  root  of  the  principal  second  moments,  all  measure 
domain  size.  It  was  then  possible  to  select  from  the  alternatives  those 
measurements  which  were  most  useful;  and  this  experience  was  used  to  reduce 
the  number  of  measurements  made  for  the  undrained  series. 

Fig. -3. 10  compares  the  Index  of  Anisotropy  calculated  from  intensity  gradient 
analysis  with  the  average  number  of  domains  in  each  image  obtained  by  con¬ 
sistency  ratio  mapping  using  eight  directions  for  the  normally-consolidated 
series  of  tests.  The  good  agreement  gave  some  confidence  in  these  new  methods 
of  image  analysis. 


3.7.3  Chord  size  analysis 

The  standard  methods  which  are  used  in  image  analysis  to  measure  features 
assume  that  all  the  features  are  comparable  in  size  and  that  most  of  them  lie 
wholly  within  the  field  of  view.  The  maps  produced  here  were  more  complicated. 
In  some  cases,  they  did  conform  to  these  assumptions;  in  other  cases,  there 
were  very  large  domains  which  overlapped  one  or  more  edges  of  the  micrograph. 
To  obtain  an  alternative  estimate  of  the  sizes  of  the  mapped  features,  the 
average  chord  size  was  calculated.  The  images  are  scanned  from  left  to  right. 
For  each  line,  the  average  chord  length  was  first  obtained  as  usual  discarding 
the  incomplete  chords  at  each  end.  Then,  if  these  incomplete  chords  were 
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.arger  than  the  average,  thev  were  Drought  hack  into  the  calculation.  Thus,  a 
global  average  for  ief t-to-r ieht  scanning  was  obtained.  The  calculation  was 
then  repeated  for  top-to-bot tom  scanning.  The  values  used  below  are  the 
averages  of  the  two  results. 


3 . S  BLACK-AND-WHITE  IMAGES 

As  an  alternative  to  intensity  gradient  analysis  and  the  mappings  based  on  it. 
black-and-white  images  were  obtained  by  thresholding  the  gray  levels  after 
Wiener  filtering.  The  threshold  was  chosen  automatically  by  comparing  each 
pair  of  adjacent  pixels  whilst  scanning  the  image  from  top  to  bottom  and  thus 
generating  a  relative  contrast  histogram  C(t).  During  the  scan,  each  pair  of 
adjacent  pixels  was  considered.  Let  /a  and  I b  be  the  gray  levels  of  such  a 
pair  of  pixels,  and  let  / b  >  fa..  Then,  for  every  gray  level  T  between  U  and 
/ b.  a  contribution  given  bv: 

dC(  T)  =  min  (T-U.  Ib-T)  (3.15). 

was  added  to  the  relative  contrast  histogram.  The  peak  of  the  relative  con¬ 
trast  histogram  based  on  all  adjacent  pairs  in  the  whole  image  was  accepted  as 
the  threshold:  pixels  brighter  than  the  threshold  were  mapped  as  particles: 
other  pixels  were  mapped  as  voids. 

Before  adopting  this  method,  several  others  were  tried.  For  simple  images,  the 
histogram  of  gray  levels  has  two  well-defined  peaks,  one  each  for  particles 
and  voids,  and  a  well-defined  and  broad  trough  between;  for  these  images,  it 
is  easy  to  select  a  satisfactory  threshold  near  the  middle  of  the  trough.  In 
our  case,  the  blurring  caused  by  the  size  of  the  point  spread  function  caused 
the  two  peaks  to  overlap  and  suppress  the  trough;  so  this  method  was 
inapplicable.  Since  the  intensity  gradient  is  expected  to  be  a  maximum  at  an 
edge,  it  was  expected  that  the  gray  level  corresponding  to  the  greatest 
intensity  gradient  would  be  acceptable  as  the  threshold;  Kohler  (1981)  reviews 
methods  based  on  this  concept.  This  concept  was  tested  using  both  3x3  and  5x5 
filters  to  calculate  mod  grad  I;  then,  for  each  gray  level,  mod  grad  I  was 
summed  and  divided  by  the  number  of  pixels  which  had  that  gray  level  to  give 
an  average  contrast  histogram;  finally,  the  gray  level  at  the  peak  of  the 
relative  contrast  histogram  was  tested  as  a  threshold.  This  method  proved  to 
be  unsatisfactory,  and,  as  noted  by  Kohler  (1981),  the  result  is  dependent 
upon  the  actual  formula  used  for  mod  grad  I.  The  idea  of  considering  only 
adjacent  pixels  was  introduced  by  Kohler  (1981);  but  he  normalized  the  con¬ 
tributions  to  the  relative  contrast  histogram  by  dividing  each  contribution 
from  each  pair  of  pixels  by  the  total  number  of  contributions  from  that  pair 
of  pixels: 


dC«c(r)  =  min  (T-I a,  Ib-T)  /  ( lb  -  U )  (3.16). 

Kohler’s  version  was  found  to  place  too  much  weight  on  the  contributions  from 
the  brighter  particles  and  was  unsatisfactory  here.  Our  decision  to  scan 
vertically  rather  than  horizontally  was  taken  after  trying  both  alternatives: 
there  was  little  difference  between  them,  and  vertical  scans  were  chosen 
because  more  particles  were  horizontal  than  vertical. 
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Fis.  3."  shows  an  lmaee  thresholded  by  this  method.  Note  the  sharpening  of  the 
image  obtained  by  Wiener  filtering  and  the  better  representation  of  the  si:e 
and  shape  of  the  voids  given  by  thresholding  after  filtering. 


3.0  POROSITY  ANALYSIS 

3.9. 1  Introduct ion 

Two  methods  were  developed  to  measure  the  local  porosity  based  respectively  on 
the  original  gray  images  and  on  those  produced  by  Wiener  thresholding. 


3.9.2  Gray  level  method 

Rv  analogy  with  X-ray  micro-analysis,  the  porosity,  p.  of  an  individual  pixel 
is  given  approximately  by: 

p  =  (  Is  -  /)/{  Is  -  /v  )  (3.17). 

where  Is  and  7v  are  the  gray  levels  corresponding  to  pure  solids  and  pure 
voids  respectively.  This  method  was  developed  for  use  in  images  in  which  the 
finest  voids  are  smaller  than  the  pixels  with  a  view  to  measuring  the  porosity 
of  each  class  of  domain.  The  results  quoted  by  Smart  and  Leng  (1990a)  for  a 
mapped  image  were  obtained  in  this  way.  If  necessary,  a  gamma  correction 
should  be  nerformed  before  making  this  analysis. 

For  the  undrained  samples,  smoothed  porosity,  p.  was  plotted  as  a  function  of 
smoothed  orientation,  T,  calculated  using  Eqn.3.13,  e.g.  Fig. 3. 8. d.  As  a  gen¬ 
eral  rule,  porosity  was  lowest  areas  aligned  in  the  direction  of  preferred 
orientation.  The  smoothed  porosity  was  also  used  for  density  field  mapping, 
see  below. 


3.9.3  Black-and-white  method 

After  Wiener  filtering  and  thresholding,  an  estimate  of  the  porosity  was  cal¬ 
culated  from  the  black-and-white  images  by  counting  the  numbers  of  pixels 
which  had  been  labelled  as  particle  and  void  respectively.  Fig. -3. 11  compares 
the  porosity  found  by  image  analysis  with  the  moisture  content  found 
gravimetrical ly;  although  not  perfect,  the  agreement  between  the  two  methods 
gives  some  confidence  in  the  accuracy  of  the  method  used  to  threshold  the 
images . 


3.10  SUBSIDIARY  ANALYSES 

3.10.1  Introduct ion 

In  addition  to  the  main  set  of  algorithms  developed  in  this  study,  several 
others  were  developed,  assessed,  or  suggested;  and  these  additional  approaches 
are  summarized  below. 
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' . 10.2  Enhanced  orientation  analysis 

Intensity  gradient  analysis  nad  originally  been  developed  in  soil  micro- 
structural  analysis  to  measure  the  anisotropy  <  f  samples  as  seen  m  scanmnu 
electron  microscopes  by  obtaining  the  polar  histogram  of  A.  These  .'.ere  often 
almost  elliptical,  whereas  Proctor  (19“’5).  whilst  working  lor  one  of  us.  had 
shown  that  some  at  least  of  the  hand  mappine  data  followed  semi—,  i rcular 
normal  distributions.  Tne  consistency  ratio  mapping  program  was  therefore 
modified  to  give  an  enhanced  orientation  analysis  by  providing  histograms 
instead  of  maps  : Leng  and  Smart.  1991).  The  polar  histograms  of  the  smoothed 
angle.  F.  give  a  stronger  and  clearer  indication  of  preferred  particle 
orientation  than  do  those  of  A ,  Fig. 3. 8. c.  These  histograms  sometimes  show  a 
rew  strong  peaks,  suggesting  a  multi-modal  semi -circular  normal  distribution: 
these  peaks  seem  to  be  in  agreement  with  subjective  examination  of  the  images. 
It  is  thought  that  the  polar  histogram  of  a  is  elliptical,  because  noise  from 
the  ends  of  the  linear  features  and  from  other  sources  has  degraded  the  'true' 
distribution.  It  is  also  possible  to  plot  the  strength  of  orientation, 
measured  for  example  by  Cor  by  the  mean  value  of  R.  against  the  direction  of 
orientation.  T.  Fig. 3. 8. b. 


3 .  10.3  Hough  transforms 

The  concept  of  multi-modal  semi-circular  normal  distributions  was  given  some 
support  by  an  alternative  analysis  using  Hough  transforms  which  was  made  for 
us  by  Luciano  da  Costa  and  Mark  B. Sandler  at  King's  College.  London  (Costa  et 
al,  1991a.  1991b). 


3.10.4  Method  of  veins  and  convex  hull 

Some  further  support  for  the  concept  of  multi-modal  semi-circular  normal  dis¬ 
tributions  was  given  by  a  novel  method  of  analysis,  called  the  method  of 
veins,  which  was  invented  by  Daisheng  Luo  (Electrical  Engineering  Department. 
Glasgow  University).  In  this  method,  the  edges  of  the  particles  are  followed 
and  recorded  as  chain  code  (Joyce  Loebl,  1985):  then,  for  each  particle  in 
turn,  the  chain  code  is  reduced  to  provide  a  measure  of  the  orientation  of 
that  particle  (Lou  et  al.,  1992).  Mr.  Lou  went  on  to  analyses  the  particles  by 
superimposing  the  convex  hull.  i.e.  the  figure  formed  by  the  enveloping 
tangents,  from  which  he  took  the  longest  diameter  as  an  indicator  of  orienta¬ 
tion.  These  two  methods  were  developed  whilst  the  top-contouring  method  was 
first  being  implemented  as  an  insurance  against  its  failure,  and  we  also  con¬ 
sidered  finding  the  principal  axes  of  each  and  every  particle;  but  all  these 
methods  would  now  seem  to  duplicate  other  findings  and  have  not  been  used  in 
the  main  analysis. 


3.10.5  Enhanced  orientation  mapping 

Rather  than  mapping  the  domains  into  predetermined  directions,  it  would  be 
preferable  to  map  them  into  the  directions  indicated  by  the  enhanced  orienta¬ 
tion  analysis.  At  the  time  of  writing,  this  is  yet  to  be  implemented:  no  dif¬ 
ficulty  is  foreseen  in  choosing  the  directions  interactively  from  the 
histogram  of  7\  but  more  experience  seems  desirable  before  attempting  to  have 
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the  directions  chosen  automatically.  it  is  thouaht  that  this  methou  .a 
enhanced  orientation  mappinc  mav  produce  an  extra  class  of  area  in  wnich  uuail 
packets  of  sub-parallel  particles  are  mixed  together  randomiv. 


3 . 1 0 . o  Semi-variogram 

Complete  semi-var ioerams  for  all  lags  up  to  one-quarter  image  size  were  calcu¬ 
lated  experimentally  for  a  few  micrographs  following  Webster  and  Oliver 
(1990).  Cross-sections  of  one  of  these  is  shown  in  Fig. -3. 9. c.  As  a  first 
approximation,  the  circular  model  semi-variogram  (Webster  and  Oliver.  1990. 
P-229)  was  fitted  to  these  results,  giving  a  range  of  18  pixels.  A  least 
squares  fit  was  used,  so  this  fit  was  objective  not  subiective.  Some  interest¬ 
ing  observations  follow. 

1.  The  range  of  the  semi-variogram.  which  represents  the  average  distance 
between  two  pixels  which  are  quite  un-related  to  each  other,  was  comparable  to 
the  size  of  the  filter  which  had  been  selected  for  top-contouring  and  con¬ 
sistency  ratio  mapping.  This  observation  increased  our  confidence  in  our 
choice  of  size  of  radius  for  these  algorithms. 

2.  The  range  of  the  semi-variogram  varied  in  direction,  thus  giving  for  the 
first  time,  a  definite  indication  of  the  ratio  of  the  axes  which  should  be 
used  if  the  circular  filter  were  to  be  replaced  by  an  elliptical  one. 

3.  Some  long  range  order  was  apparent.  This  observation  led  to  the  develop¬ 
ment  of  the  analysis  described  below  under  Density  Field  Mapping. 

4.  The  nugget  variance,  i.e.  the  extrapolated  variance  at  zero  lag,  was  zero, 
suggesting  that  the  magnification  could  have  been  lowered  without  serious  loss 
of  fine  detail  and  with  a  gain  in  the  size  of  the  area  examined. 

These  observations,  which  were  made  recently  by  Dr.Leng,  require  further 
development  and  are  the  subject  of  a  research  application  by  her. 


3.10.7  Density  field  mappings 

After  discovering  long  range  order  in  the  semi-variograms,  the  variability  of 
the  density  field  of  the  undrained  samples  was  analyzed.  Instead  of  density 
itself,  the  local  porosity  was  calculated  by  the  gray  level  method  within  the 
same  large  uniform  filter  as  was  used  for  consistency  ratio  mapping.  An 
intensity  gradient  analysis  was  made  of  the  porosity  field  and  compared  with 
the  results  from  the  original  intensity  gradient  method  and  from  enhanced 
orientation  analysis.  To  obtain  the  scale  of  the  variation  of  porosity,  the 
porosity  field  was  segmented  at  the  mean  value  for  the  sample,  and  a  chord 
size  analysis  made. 


3.10.8  Further  porosity  analysis 

After  images  have  been  Wiener  filtered  and  thresholded,  the  voids  are  clearly 
distinguished  from  the  particles.  Two  methods  have  been  considered  to  extract 
further  information  about  the  voids.  In  one  method,  the  voids,  which  are  long 
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and  tortuous,  could  be  sub-divided  at  their  'necks'  and  the  pieces  analyzed 
.ndividuai ly.  Alternatively,  ihe  particles  could  be  dilated  repeatedly  by 
ncreasing  amounts,  thus  giving  a  pore  size  distr lbut ion. 
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Appendix  3.1  LARGE  UN I FORM  FILTERS 


-'he  i  liters  useu  here  are  necessarily  !.  ir.se .  typical  iy  2?  pixels  radius,  con¬ 
taining  about  1300  pixels.  However,  the  algorithms  used  .ire  additive.  There¬ 
fore.  uniform  filters  were  specified  to  limit  the  work  required  as  the  filter 
traverses  the  image.  At  each  pixel  in  a  line,  the  trailing  edge  of  the  filter 
is  dropped  and  a  new  leading  edge  is  added.  Thus,  for  example,  only  i.2  pixels 
t round  the  perimeter  ire  involved,  a  reduction  of  more  than  *<o  The  shape  of 
the  filter  is  defined  by  the  half  chord  lengths,  e.g  RHS(l-) .  .here  V  is  the 
distance  of  the  chord  from  the  (.enter  of  the  filter,  -'ircular  filters  are  used 
to  minimize  spurious  orientation  effects.  The  filter  radius  is  specified  in 
pixels,  but  0.4  is  added  because  this  was  found  to  improve  the  circularity  of 
-ircies  up  to  ~  pixels  radius.  For  a  nominal  radius  of  i.  the  eight  nearest 
neighbors  are  used:  arid  for  a  nominal  radius  of  0.  the  four  nearest  neighbors 
are  used.  In  the  Fortran  versions,  the  two  sides  of  the  filter  are  defined 
separate  1 y : 


LHS(-P)  =  RHS(P) 


(3.  IS)  . 


.0  that  shapes  such  as  inclined  ellipses  could  be  used  if  necessary.  Square, 
rectangular,  and  straight-line  filters  are  also  possible. 


For  moving  from  line  to  line,  the  top  edge  is  dropped  and  a  newi  bottom  edge  is 
added.  The  top  and  bottom  of  the  filter  are  defined  by  TOP(L’)  and  BOT(U). 
respectively,  where  l’  is  the  distance  of  the  chord  from  the  center.  A  separate 
set  of  registers  is  maintained  for  the  first  pixel  in  each  line,  and  it  is 
these  which  are  updated:  then  their  contents  are  copied  to  the  registers  for 
the  filter  which  traverses  the  line. 


Three  schemes  were  considered  to  avoid  unnecessary  tests  within  the  double 
loop  needed  to  process  the  whole  image. 

a.  In  the  first  version  of  the  code:  the  image  array  exactly  equalled  its 
nominal  size:  and  the  first  filter  on  each  line  was  set  up  without  its  leading 
edge.  Then  the  operations  within  the  loop  were:  add  leading  edge:  calculate: 
drop  trailing  edge.  A  similar  procedure  was  used  for  the  'leading  column’  fil¬ 
ter. 

b.  In  the  present  code:  the  image  array  is  extended  by  one  line  and  by  one 
column:  (the  extra  pixels  are  set  to  zero  for  definiteness):  and  the  first 
filter  on  each  line  is  set  up  naturally.  Then  the  operations  within  the  loop 
are:  calculate;  drop  trailing  edge;  add  leading  edge.  At  the  last  pixel  in 
each  line,  adding  the  leading  edge  is  unnecessary;  but  it  takes  less  time  than 
would  repeated  testing  to  see  whether  the  end  had  been  reached,  and  it  is  the 
simplest  method  to  code. 

c.  If  the  image  array  cannot  be  extended,  then  the  loop  could  be  shortened  by 
one  pixel,  and  the  calculations  for  the  last  pixel  could  be  made  outside  the 
loop.  This  seems  to  demand  the  use  of  INCLUDE  statements  for  the  repeated  por¬ 
tions  of  code  but  would  perhaps  make  the  development  of  new  algorithms  easier. 

When  using  large  filters,  it  is  important  to  avoid  leaving  large  unprocessed 
borders.  Smart  and  Leng  (1990a)  showed  how  to  do  this  by  reflection:  and  Smart 
et  al.  (1991)  have  a  more  general  discussion. 
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Appendix  3.2  MEASUREMENTS  FROM  DOMAIN  MAPS 

The  following  quantities  were  measured  for  all  the  domains  and  random  clusters 
•thich  had  been  mapped: 


! 

XREF 

Domain  reference  point,  bottom  right-most  pixel 

T 

VREF  ; 

3 

ID 

Identification  number 

4 

PARENT 

Identification  number  of  parent  region  surrounding  domaii 

5 

HOLES 

Number  of  holes  in  domain 

BACKGROUND 

Flag  indicating  domain  is  a  hole  inside  another  domain 

CONTACT 

Flag  indicating  domain  touches  (overlaps)  edge  of  image 

0 

XMIN  ) 

Limits  of  domain  in  image  coordinates 

[) 

XMAX  ) 

:o 

YMIN  ) 

1 1 

VMAX  ; 

12 

HFERET 

Feret  diameter,  horizontal 

13 

VTERET 

vert ical 

14 

AFERET 

at  45  degrees 

]  5 

BFERET 

at  135  degrees 

16 

HPROJ 

Total  projected  edge  length,  horizontal 

r 

VPROJ 

vert ical 

18 

PERIMETER 

Perimeter 

19 

.AREA 

Area 

20 

XCEN  ) 

Centroid 

21 

YCEN  ) 

n 

MM  IN  ) 

Principal  second  moments  of  area 

23 

MMAX  ) 

24 

ANGLE 

Orientation  of  principal  axis  of  domain 

25 

CIRCULARITY 

4  *  P I  * AREA/PER I METER*  *  2 

The  following  quantities  were  calculated  subsequently: 

1  Mean  Feret  diameter,  estimates  domain  size 

2  Equivalent  circular  radius,  estimates  domain  size 

3  SQRT(MMIN*MMAX) ,  estimates  domain  size 

4  Minimum  Feret  diameter/Maximum  Feret  diameter,  estimates  domain  shape 

5  MMIN/MMAX.  estimates  domain  shape 

6  (HPROJ+VPROJ )/( 2*mean  Feret  diameter),  estimates  indentation 

n  ( HPROJ+VPROJ )/(MM I N+MMAX) ,  estimates  indentation 
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Fig.  3.1  Typical  scanning  electron  micrograph. 

Fig.  3.1. a  Original  image. 

56  um  square.  Natural  orientation:  up-down  =  vertical: 
across  =  horizontal. 

(Fig.  3.1.b  follows  on  the  next  page.) 


I 


Fig.  3.2  Simplified  flow  chart,  Glasgow 


1)  Image  Acquisition 

2 )  Co mpute  Courier  Transform 


2)  Search  for  Interference  Peaks 

and  Mask 


4)  Determine  Signal  to  Noise 


5) 


Image  Reconstruction 
Wiener  filtering 


Compute  relative 
contrast  histogram 


7) 

8) 

9) 

10) 


_ Threshold  image 

Porosity  Computations 
Domain  Segmentation 

Porosity  Computation 
for  each  domain  class 


Fig.  3.3  Simplified  flow  chart,  UEA 
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Fig.  3.4  Image  analysis  facility  in  UEA. 

A.  PCs  with  Synapse-Sprynt  Frame  Grabber  and  Semper. 

B.  PC  with  Synergy  Slow  Scan  Frame  Grabber  and  Semper. 

C.  Sun  Workstation  with  additional  hard  disk. 

Exobyte  tape  streamer  and  Semper. 

D.  PC  serving  laser  printer  also  connecting  Sun  to  netwoik. 

E.  PC  with  optical  disk  archive  and  tape  streamer. 

F.  PC  holding  master  progran.s  and  connecting  Sun  to  network. 

G.  Interface  between  S800  SEM  and  Synergy  Frame  Grabber 

H.  Control  switch  to  select  various  other  image  types. 

LINK.  Link  Analytical  digitiser  for  S450  SEM. 
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w.  Windows  in  Sun  Workstation 
MO  Master  transputer 
W.  Worker  transputers 
GO  Graphics  transputer 
S  Image  screen 
D  Disk  store 
T8  Independent  transputer 
Links  are  two-way,  except  to  S 


Fig.  3.5  Meiko  multi-processor  system  at  Glasgow. 
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S1.S2  Image  screens;  M  PC  monitor;  Cam’  Camera: 

G300  Transputer,  graphics  chip,  2  MB  work’g,  2  MB  video: 
86/7  CPU,  2  MB  =  work’g,  system,  cache,  d:  ram  disk,  etc.; 
Semp  Semper,  framestore; 

A,B  Floppy  drives;  Keyb  Keyboard;  C  hard  disk. 


Fig.  3.6  Hybrid  dual -processor  system  at  Glasgow.  Nimbus. 


Fig.  3.7  Example  of  Wiener  filtering. 

Fig.  3. 7. a  Original  image. 

Fig.  3.7.b  Thresholded  after  Wiener  filtering. 


Fig.  3.8  Histograms  (graphs)  vs.  Orientation. 

Fig.  3. 8. a  Amount  of  orientation.  A,  intensity  gradient  analysis. 
Fig.  3.8.b  Strength  of  orientation,  average  mean  vector. 

Fig.  3.8.c  Amount  of  orientation,  T,  consistency  ratio  mapping. 
Fig.  3.8.d  Average  porosity. 
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40  60  SO  100 

RADIUS  (pixels) 


Fig.  3.9  Choice  of  radius  of  filter. 

Fig.  3. 9. a  Modified  chi-square,  for  Fig.l. 

Fig.  3.9.b  Consistency  ratio,  for  a  low  magnification  imag 
Fig.  3.9.c  Semi-variogram,  circular  model  fitted,  for  Fig. 
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Number  of  domains  per  image 

Fig.  3.10  Index  of  Anisotropy  vs.  Number  of  domains  (SEM  nc). 
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Image  Analysis  Porosity 


Fig.  3.11  Image  analysis  porosity  vs.  Moisture  content  (SEM  nc). 

The  gravimetric  moisture  content  is  thought  to  be  wrong  for  4  samples. 
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4.1  INTRODUCTION 

Mitchell  was  the  first  to  apply  the  technique  of  optical  microscopy  to  the 
investigation  of  soil  structure  in  geotechnical  engineering  in  1956  (Mitchell. 
1956).  Much  use  has  been  made  of  optical  microscopy  since  then  (e.g.  Wu.  1958, 
Smart.  1966a:  Tchalenko.  1967).  The  usual  technique  and  the  one  which  was  used 
here  employs  thin  sections.  These  were  prepared  conventionally  from 
impregnated  'half -samples*  similar  to  those  used  for  scanning  electron  micro¬ 
scopy  (see  previous  Chapter). 

The  technique  for  observation  of  micro-structural  features  in  thin  sections  of 
clay  is  based  on  the  optical  properties  of  clay  particles.  A  platy  kaolinite 
particle  is  optically  anisotropic  and  possesses  three  approximately  mutually 
perpendicular  optical  axes:  two  almost  in  the  plane  of  the  particle,  which  are 
called  section  axes  here,  and  one  almost  at  right  angles  to  this  plane.  The 
refractive  indices  in  the  directions  of  the  section  axes  are  approximately 
equal  and  significantly  different  from  that  in  the  third  axial  direction.  When 
a  single  particle  is  viewed  between  crossed  poiars,  it  may  or  may  not  exhibit 
extinction  and  illumination;  but  when  a  thin  section  is  examined,  it  transmits 
a  variable  intensity  of  light  depending  on:  the  orientation  of  the  optical 
axes  of  the  particles  in  the  thin  section,  the  refractive  indices  of  the 
mineral,  the  thickness  of  the  thin  section,  the  wave  length  and  the  intensity 
of  the  light.  This  transmitted  intensity  varies  as  the  thin  section  is  rotated 
with  respect  to  the  direction  of  polarization.  For  a  certain  thin  section, 
assuming  that  its  thickness  is  uniform,  under  a  certain  light,  the  intensity 
of  light  transmitted  through  the  thin  section  depends  only  on  the  orientation 
of  the  optical  axes  of  the  particles  in  the  thin  section  and  the  related  posi¬ 
tions  of  the  optical  axes  of  the  particles  and  of  the  direction  of  the  light. 
Thus,  if  kaolinite  particles  are  arranged  in  a  random  style  within  a  sample,  a 
thin  section  in  any  direction  through  the  sample  does  not  exhibit  extinction 
or  illumination,  since  the  optical  axes  of  the  particles  are  also  distributed 
randomly.  When  such  a  thin  section  is  rotated  between  crossed  poiars,  the 
intensity  of  transmitted  light  remains  constant  throughout  the  360*.  If  the 
particles  are  oriented  with  their  planar  surfaces  parallel,  extinction  and 
illumination  will  occur  in  thin  sections  cut  from  the  sample  in  particular 
directions.  When  a  thin  section  is  cut  parallel  to  the  orientation  planes  and 
viewed  between  crossed  poiars,  no  change  shows  in  illumination  with  rotation, 
since  the  refractive  indices  in  the  directions  of  the  section  axes  of  the  par¬ 
ticles  are  approximately  equal;  but,  when  a  thin  section  is  cut  normally  to 
the  plane  of  preferred  orientation,  four  stages  of  extinction  and  illumination 
occur  between  crossed  poiars,  since  the  optical  axes  of  the  particles  lie  in 
preferred  directions.  Hence,  measurement  of  the  transmitted  light  intensities 
may  be  used  as  an  indirect  measurement  of  the  particle  orientation.  This  Chap¬ 
ter  gives  a  review  of  a  previous  method  and  introduces  three  new  and  improved 
methods. 
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-.2  POLARIZING  MICRO-PHOTOMETRIC  ANALYSIS.  PMPA 

Previous  work  has  used  Polarizing  Micro-Photometric  Analysis.  PMPA.  ;  or  .. 
microscope  with  crossed  polars  and  a  thin  section  between  them,  for  the 
intensity  of  light,  the  transfer  function  is  (e.g.  Smart.  1966a;: 

1  =  0.5  sin2  10  sin2  o/2  (4.1). 

where  0  is  the  ingle  between  the  section  axes  and  the  axes  of  the  poiars:  J  is 
the  phase  difference  after  the  light  r.as  passed  through  a  thin  section,  which 
is  a  function  of  the  thickness  of  the  thin  section,  d.  and  the  wave  leneth  of 
the  light.  1.  Therefore,  for  areas  not  dark  in  ordinary  light,  the  diagnosis 

i  s : 

-  dark  for  all  0 : 

r  for  all  wavelengths: 

r  when  section  tilted:  Isotropic. 

!  otherwise:  Optic  axis  x  section. 

l 

*-  only  for  wavelength  =  1:  d  =  kl/Sn .  k  -  1.2 . 

L  4  extinction  positions:  Extinction  directions  -1-  section  axes. 

(per  revolution) 


As  mentioned  above,  it  may  be  found  that  the  greater  the  difference  between 
illumination  and  extinction,  the  greater  the  degree  of  the  orientation  of  par¬ 
ticles.  So  with  a  microphotometer,  quantitative  results  are  obtained  by 
measuring  the  intensities  of  the  polarized  light  transmitted  through  a  section 
in  the  illumination  and  extinction  positions.  The  ratio  of  the  two  measure¬ 
ments.  or  some  other  combination  of  factors  dependent  on  them,  is  an  indica¬ 
tion  of  the  degree  of  preferred  orientation. 

Smart  (1966a)  defined  an  Anisotropy  Index.  Mi .  as: 


Ml 


Im  ax  ~  Iain 
I  max  +  la  i  n 


(4.2), 


where  /  is  the  intensity  of  the  light  transmitted.  In  theory,  the  value  of  the 
anisotropy  index  Mi  ranges  from  0.  which  indicates  a  random  arrangement  of 
particles,  to  1,  which  indicates  an  ideal  orientation  arrangement  of  parti¬ 
cles.  In  practice,  because  of  slight  asymmetry  of  the  measurements,  an  average 
of  four  maxima  and  an  average  of  four  minima  are  used.  The  difference  between 
the  formula  for  Anisotropy  Index,  Mi,  here  and  Index  of  Anisotropy,  /a,  used 
for  scanning  electron  microscopy  was  deliberate,  as  it  was  thought  that  the 
physical  conditions  of  measurement  were  different. 


By  means  of  the  microphotometer,  the  average  orientation  direction  of  the  par¬ 
ticles  can  be  determined  too.  A  full-wave  compensator  could  be  used  to  locate 
the  direction  of  preferred  orientation  by  comparing  the  colors  seen  with  the 
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.olors  .f  failure  planes.  However,  '.here  :  i  mall  risk  of  an  error  ol  no’*  in 
any  male  observation. 

In  the  observations  described  here,  the  objective  was  xlO  and  the  field  of 
"iew  was  1.5  mm  diameter.  In  general,  live  fie.ds  of  view  were  observed  from 
one  thin  section  from  each  triaxial  test  sample. 


4.3  POLARIZING  MICRO-PHOTOMETRIC  MAPPING  I,  POLMAP 

■'•hen  viewed  between  crossed  polars.  the  matrix  of  clay  soil  often  shows  dif¬ 
fuse  light  and  dark  patches  which  represent  information  which  has  been  diffi¬ 
cult  to  extract  by  conventional  microscopy.  By  using  a  special  type  of  optical 
microscopy  and  computer  analysis  system,  the  following  new  method  has  been 
developed. 

It  is  known  that  extinction  occurs  only  when  the  section  axes  are  parallel  to 
or  perpendicular  to  the  axes  of  polars  when  a  thin  section  is  viewed  between  a 
crossed  polars;  and,  when  a  quarter-wave  plate  is  used  in  such  a  position  that 
the  angle  between  the  axes  of  the  quarter-wave  plate  and  the  axis  of  the 
polarizer  is  45°.  the  transfer  function  becomes: 

I  =  0.5  ( 1  -  sin  20  sin  J)  (4.3). 

In  this  case,  full  extinction  occurs  when  the  angle  between  the  section  axes 
and  the  axi.,  of  the  polarizer  becomes  45°  and  6  is  sin_1  +  l  or  this  angle  is 
135°  and  <5  is  sin"l 2-l. 

Most  of  the  measurements  discussed  here  were  made  using  a  Leitz  Dialux-pol 
microscope  in  which  the  polarizer,  the  analyzer,  and  a  slot  for  a  compensator, 
are  all  linked  together  by  an  arm  and  can  be  rotated  synchronously  (see  Smart 
and  Tovey  19S2  Fig. 2.1).  The  Leitz  microscope  used  here  has  another  slot  which 
is  fixed  in  the  body  of  the  microscope,  and  which  must  not  be  used  here.  Some 
measurements  were  made  using  a  Prior-Swift  MPB124  Polar  Microscope  in  which 
similar  synchronous  rotation  was  achieved  by  internal  gearing. 

To  explain  the  method,  Fig. 4.1  shows  a  failure  plane  running  diagonally 
through  a  field  most  of  which  is  oriented  up-and-down  the  page;  the  field  of 
view  is  0.47  mm  square. 

In  Fig. 4.1. a.  the  polars  were  crossed  "on  the  square’,  the  synchronizer  angle 
being  90°.  so  the  failure  plane  is  bright;  and  dark  areas  are  areas  in  which 
the  clay  particles  are  aligned  either  horizontally  or  vertically  (with  respect 
to  the  field  of  view).  This  micrograph  was  digitized  and  denoted  by  x2.  A  cor¬ 
responding  micrograph,  xl,  was  digitized  with  the  polars  crossed  ’on  the 
diagonal’,  the  synchronizer  angle  being  45°;  in  this  micrograph  the  contrast 
was  reversed.  In  the  first  stage  of  analysis,  these  two  microphotographs  were 
compared  to  distinguish  between  areas  oriented  ’on  the  square’  and  ’on  the 
diagonal’.  Each  pixel  was  coded  as  follows: 

1  if  brighter  in  xl  than  x2,  on  the  square; 

2  if  brighter  in  x2  than  xl.  on  the  diagonal; 

0  otherwise,  i.e.  undecided. 


M'ter  this  preliminary  coding.  four  more  digitized  micrographs  ‘ere  used  to 
distinguish  horizontal  from  vertical  ind  down- i eft  i  rom  down- 


mi  crop  rapns  cere  taken  with  a  quarter- wave  plate 
slot  and  rotating  synchronous ly  with  the  poiars. 


Inserted  in 


right.  All  these 
the  compensator 


in  the  second  stage  of  analysis,  only  pixels  _eded  2  were  considered,  anu  two 
micrographs  taken  with  the  synchronizer  set  at  ,J0°  and  0”  were  used.  These  two 
micrographs  are  denoted  q2  and  q4  and  are  shown  in  Figs. 4.  i.e  and  d.  respec¬ 
tively.  In  one  of  these  micrographs,  the  brighter  areas  are  down-left,  and  in 
the  other  they  are  down-right:  a  pragmatic  method  of  deciding  which  is  which 
is  to  experiment  with  failure  planes  in  which  the  particle  alignment  may  be 
assumed  safely.  Thus,  the  pixels  which  had  originally  been  coded  2  were 
"ecoded  as  follows: 

2  if  brighter  in  q2  than  q4,  i.e.  down-left: 

4  if  brighter  in  q4  than  q2.  i.e.  down-right: 

0  otherwise,  i.e.  undecided. 

Similarly,  using  micrographs  ql  and  q3.  with  the  synchronizer  -et  at  135°  and 
45°.  respectively,  the  mapping  can  be  completed.  The  result  is  shown  in 
Fig.4.t.b.  where  the  two  diagonal  classes  have  been  ruled  and  the  up  and 
across  areas  have  been  left  blank. 

During  analysis,  a  convenient  method  of  presenting  the  results  is  10  color  the 
pixels  as  fol  low's: 


1 

1 

3 

4 

0 


horizontal , 
down- left , 
vert ical , 
down-right , 
undecided. 


red: 

yellow: 

green; 

blue; 

black. 


The  colors  were  chosen  because  thin  sections  such  as  the  one  used  here  fre¬ 
quently  show  yellow  and  blue  for  the  pixels  coded  2  and  4  respectively  when 
viewed  in  white  light  with  the  synchronizer  set  at  90°  and  using  a  full-wave 
compensator.  Under  those  conditions,  the  other  areas  tend  to  be  pink,  leading 
to  green  vertical  and  red  horizontal.  In  theory,  for  the  analytical  method 
used  here,  monochromatic  light  with  a  wave-length  matched  to  that  of  the 
quarter  wave  plate  should  be  used:  but  in  practice,  white  light  was  found  to 
be  more  convenient  and  to  give  similar  results. 

The  overall  anisotropy  of  the  field  of  view  may  be  assessed  by  considering  the 
number  of  pixels  mapped  into  each  direction:  undecided,  horizontal,  down-left, 
vertical,  down-right.  Following  the  PMPA  method  above,  an  Anisotropy  Index. 

.4g,  obtained  by  this  method,  polarizing  micro-photometric  mapping  (PMPM),  may 
be  defined  as  follows: 


.-12 


(H  +  n  -  (L  +  R)  | 
H  +  V  +  L  +  R 


(4.4) 


where  H .  T.  L  and  R  are  respectively  the  number  of  pixels  mapped  into  the 
directions:  horizontal,  vertical,  down-left  and  down-right.  The  Anisotropy 
Index.  .42.  obtained  from  POLMAP  can  be  compared  with  that  given  by  polarizing 
micro-photometric  analysis  (PMPA),  ,4i .  Usually,  for  the  same  area,  .4i  was  a 
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.  i  1 1 1  i .r.ver  '.ban  .u.  ;  hese  two  values  of  .1  are  r.ot  trict  i  v  .amparable. 

>e cause*  the  t'MPA  measurement  has  !ven  t.a.en  over  a  ider  'held  ;:t  •.  lew  than 
the  FOI.MAP  measurement.  Further .  the  direct  ions  usees  far  ■OLMAP  vere  zhosen 
irbit  rarily  beforehand.  ..hereas  the  direct  ions  used  for  ."IP  A  acre  the  actual 
observed  direct  sons  of  maximum  and  minimum  preferred  orientation. 

rh.e  sizes  of  the  patches  mapped  by  l 'OLMAP  .vere  measured  usinu  the  ..hord  sire 
algorithm  nich  tvas  described  in  the  previous  •  hapter. 

in  general.  Four  fields  from  one  thin  section  vere  measured  for  each  triaxiai 
ample  usinu  F’OLMAP:  all  the  measurements  quoted  lie  low  .vere  made  on  fields 
' .  4”  mm  square. 

Three  torsions  of  another  new  method  are  indicated  below. 


4.4  POLARIZING  MICRO- PHOTOMETRIC  MAPPING  II 
4.4.1  Introduct ion 

After  establishing  the  POLMAP  algorithm,  a  set  of  algorithms  which  obtain  the 
orientation  of  each  pixel  was  developed.  Of  the  versions  described  below.  TOL- 
MAP3  is  the  most  economical:  but  the  digitized  micrographs  required  for  POL- 
MAP4  had  already  been  acquired,  so  this  version  was  used  for  the  main  analy¬ 
sis. 


4.4.2  Po 1 map4 

When  a  parallel  beam  of  monochromatic  light  passes  through  a  polarizer,  a  thin 
section,  a  quarter-wave  plate,  and  an  analyzer  which  is  crossed  with  the 
polarizer,  as  mentioned  above.  I,  the  transfer  function  of  the  intensity  of 
the  light  is  as  follows: 

I  =  0.5  11  -  sin  20  sin  J)  (4.5). 

where  0  now  refers  to  the  inclination  of  one  of  the  section  axes  of  a  pixel. 
Then,  suppose  that  the  synchronizer  is  set  in  turn  at  0°.  45°,  90°.  and  135°. 
then,  the  intensity  of  the  light  at  each  setting  may  be  calculated  from  the 
following  functions. 


Iq 

=  0.5 

(  1 

-  sin 

20  sin 

S) 

cr 

m 

(4.6). 

Us 

=  0.5 

(  1 

-  cos 

20  sin 

5) 

=  q3 

(4.7). 

1 90 

=  0.5 

(  1 

+  sin 

20  sin 

3) 

=  q4 

(4.3). 

/l  3  5 

=  0.5 

(  1 

+  cos 

20  sin 

6) 

=  ql 

(4.9), 

where  ql,  q2,  q3  and  q4,  denote  the  corresponding  micrographs.  Therefore: 

bo  -  Jo 


tan  20 


1 13S  -  Us 


(4.10), 
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'■  i  n  j  -  .  j  o  -  _  \j  )  *■  *  -t ...  • 

Chus .  the  direction  f  preferred  crxentat  son  .;t  .•  ich  pi:--;  -an  K-  .’-tamed. 
..vcept  that  ;»  was  labelled  undee  iued  :  f  :  he  numerator  tna  denominator  ot  tan 
1J  were  both  snail,  rim.  new  method  a  as  called  TiLMAP4. 

"ig.4.2  snows  a  poiar  histogram  for  one  field  of  : e’.v  untamed  u.>ir.g  this 
.net  hod. 

for  mapping,  the  degrees  .’.ere  divided  into  groups,  men  acre  horizontal 
22.5"  and  I5~. 5"-l80“  '  .  down-left  {  22. 5“ -6~ .  5 "  ) .  vertical  l  c“ .  5°- i  i  2 .  ?° ! .  and 
Iown-nght  ;  !  12.5"-15“.5“ ) .  Comparison  with  POLMAP  gave  satisfactory  result.;. 
If  desired,  the  degrees  could  he  divided  into  a  larger  number  of  groups,  -say 
or  12.  for  a  finer  mapping. 

Taking  the  assumed  intensity  at  each  pixel  as  a  unit  vector,  the  Consistency 
Patio.  C.  was  calculated  from: 


A  = 

2  cos 

20 

(4.12) 

y  - 

1  sin 

20 

(4.13) 

R  = 

>  ( A- 

+  5"-  ! 

(4.14) 

C  = 

R/S 

(4.15) 

T  = 

i  atan 

i  y/x) 

(4.16) 

where  S  is  the  number  of  decided  pixels,  f  was  labelled  undecided  if  .V  and  V 
were  both  small.  The  values  of  Consistency  Ratio  obtained  were  almost  always 
higher  than  the  corresponding  values  of  Anisotropy  Index. 

For  the  normally-consolidated  undrained  series  of  triaxial  tests.  POLMAP4  was 
used  in  addition  to  POLMAP. 


4.4.3  Pol map 3 

P0LMAP3  comes  from  the  same  theory  as  P0LMAP4  and  produces  the  same  results 
from  only  three  micrographs,  when  the  synchronizer  is  set  in  turn  at  0*.  120°. 
240°.  the  intensities  of  the  light  through  the  thin  section  are  determined  by 
the  following  functions  respectively: 


io 

=  0.5  (1  - 

sin  20  sind) 

(4.17). 

1 120 

=  0.25  (2  + 

(sin  20+  *3  cos 

20)  sind) 

(4.18). 

id  4  0 

=  0.25  (2  + 

(sin  20  -  y"3  cos 

20)  sind) 

(4.19). 

Therefore: 

tan  20  = 

ii  2  0  +  I2  4  0  -  2  io 

(4.20), 

J 3  (  il  20  -  i240  ) 

sin  6  =  2/3  /{(ii2o  +  Ii 40  - 

■  2  io  ) 2  +  3  (  i  120 

-  i2  4  0  ) 2  } 

(4.21). 
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i  1  >2 . 4 .  :  a i so  ■'nows  •. he  polar  histogram  ;  •. he  -arae  :u‘id  of  i e»  obtained 

■Jsma  P0LM.XP3 .  .vnica  has  •  ery  rood  agreement  .  1 1  a  .ne  result  tVorn  :’'L‘LAP4. 

P.v  the  same  method,  '.he  -insistency  Ratio  muv  uso  ne  calculated  from  POLMAP 3 . 

P0LMAP3  requires  three  micrographs,  qo  .  q;  :o  .  md  q:4o  lor  the  :-r  ior-Sw  i  ft 
MPB124  Polar  Microscope.  For  the  leitr  Microscope-  three  micrographs  qo .  q • ; j . 
and  quo.  may  he  used,  because  q:-u;  is  unavailable  md  ..a  =  4  o . 

.'.It hough  POLMAP3  requires  fewer  micrographs  than  either  POLMAP  or  1  •'.’LM\P4 .  the 
data  required  for  10I.MAP4  had  already  been  measured  for  POLMAP.  rtvas.  it  was 
convenient  md  saved  time  to  use  POLMAP4  here. 


4.4.4 


Po 1 mapb 


Since  lu o 
have  been 
field,  so 


=  /:4o.  etc..  POLMAP 3  generalizes  readily  to  PCLMAPb.  which  would 
more  accurate.  But  P0LMAP6  would  have  required  six  micrographs  per 
it  was  not  implemented. 


ooooooooooo 
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f  ig.  4.2  ijolar  Histograms  for  one  field  0.47  mm  square. 

Full  line:  rrom  P0LMAP4.  Broken  line:  from  P0LMAP3 .  Oriented  naturally:  up 
down  -  vertical:  across  =  horizontal. 


chapter  5 


NORMALLY -CONSOL  I  DATED  DRAINED  TESTS 


5.1  INTRODUCTION 

The  new  methods  of  micro-analysis  were  first  applied  to  a  series  of  41 
normally  consolidated  drained  triaxial  tests  of  Speswhite  kaolin  after 
preparation  in  the  medium  pressure  filter  presses.  In  brief:  the  slurry  was 
mixed  at  140  %  moisture  content,  i.e.  303  T  Liquidity  Index:  consolidated 
uniaxial  ly  to  700  kPa:  consolidated  triaxial  lv  to  "’00  kPa:  strained  at  5  %  per 
day:  impregnated  with  Araldite  AY18;  and  prepared  for  microscopy. 


5 • 3  STRESS-STRAIN  BEHAVIOR 

Fig.  5.1  illustrates  the  deviator  stress  vs.  axial  strain  for  typical  samples 
from  the  normally  consolidated  undrained  series,  the  results  being  much  as 
expected,  i.e.  a  steady  rise  of  deviator  stress  to  a  peak  value  at  say  16  Z 
axial  strain,  at  about  which  strain  failure  occurred,  there  being  a  small 
amount  of  variability  from  sample  to  sample.  Fig.  5.2  shows  the  volumetric 
strain  vs.  axial  strain  for  the  whole  normally  consolidated  series,  each  point 
being  the  end  point  of  one  test,  except  at  the  origin,  where  there  were 
several  control  samples.  As  expected,  the  samples  consolidated  throughout  the 
tests.  The  rate  of  consolidation  was  faster  at  first  and  then  slowed  con¬ 
siderably.  there  being  a  fairly  strongly  marked  change  of  behavior  at  about  8 
”  axial  strain.  Also,  somewhat  unexpectedly  perhaps,  the  points  during  the 
early  stages  of  the  tests  were  not  far  removed  from  the  Ko-line,  i.e.  the 
lateral  strain  was  almost  negligible  up  to  about  3  %  axial  strain  (e.g. 

€i  ■  ~  - 1 . 5  ~  when  £a-  =  -S-%,  but  £i  ~  -7-°5  when  6a  ■  =  •  20  •  %  ). 


5.3  VISUAL  EXAMINATION 

'.Vhen  viewed  between  crossed  polars  with  or  without  a  full-wave  plate,  the 
kaolin  often  showed  diffuse  light  and  dark  patches,  although  some  fields  were 
more  uniform  or  had  more  obvious  horizontal  preferred  orientation  than  others. 

Visual  examination  of  the  scanning  electron  micrographs  showed  that  the  struc¬ 
ture  of  the  samples  was  composed  mainly  of  diffuse  domains,  which  tended  to 
merge  together.  As  indicated  below,  there  were  a  few  random  clusters;  but 
there  were  virtually  no  inter-domain  voids  in  these  samples.  Thus,  the 
hypothesis  that  shear  deformation  is  accommodated  by  slip  of  domains  into  the 
inter-domain  voids  is  irrelevant  here;  but  there  could  still  have  been  slip  of 
the  domains  into  the  looser  and  presumably  weaker  random  areas. 

Failure  planes,  when  they  occurred,  were  slickensided  as  usual  and  cut  through 
the  previously-existing  structure  non-conformably;  however,  attention  here  is 
concentrated  on  the  ’pre-peak’  deformation. 
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rhe  measurements  discussed  below  were  taken  on  !  I  elds  irl  lows: 

FMPA  1.5  mm  d  i  ame  t  e  r 

POLMAP  0.4~  mm  square 

5EM  0.025  ir,m  square 


5.4  PREFERRED  ORIENTATION 

lie.  5.3  shows  the  preferred  orientation  vs.  axial  strain  obtained  from  the 
polarizing  micro-photometr ic  analysis,  PMPA.  Here  and  in  the  graphs  following, 
each  point  is  the  end  point  of  one  test.  The  results  in  this  figure  are  the 
moduli  of  the  algebraical  average  of  the  five  values,  of  orientation  obtained 
from  the  five  fields  measured  by  this  technique,  and  they  are  expressed  as 
inclinations  from  the  horizontal.  Fig.  5.4  shows  the  corresponding  values 
obtained  from  scanning  electron  microscopy  using  intensity  gradient  analysis. 
\1 1  the  values  in  both  graphs  are  low:  i.e.  the  preferred  orientation  started 
close  to  the  horizontal  and  remained  close  to  the  horizontal  throughout.  Thus, 
the  hypothesis  that  the  particles  would  gradually  rotate  into  the  direction  of 
the  failure  plane  is  eliminated. 


5.5  STRENGTH  OF  ANISOTROPY 

Figs.  5.5,  5.6,  and  5.'7  show  the  strength  of  anisotropy  vs.  axial  strain 
measured  by  PMPA,  POLMAP.  and  SEM,  respectively.  Each  point  is  the  arithmeti¬ 
cal  average  of  Anisotropy  Index  for  5.  4,  and  24  fields,  respectively. 

In  all  three  graphs,  the  .Anisotropy  Index  starts  at  a  moderate  value  and  ends 
at  a  moderate  value.  This  is  in  accordance  with  the  hypothesis  that,  during 
shear  deformation,  strongly  anisotropic  samples  become  less  anisotropic  and 
weakly  anisotropic  samples  become  more  anisotropic;  these  samples  appear  to  be 
on  the  borderline. 

All  three  graphs  appear  to  show  a  further  effect.  To  illustrate  this,  quad¬ 
ratic  regressions  have  been  added  to  the  first  two  figures,  and  linear  regres¬ 
sions  have  been  added  to  the  early  and  late  stages  of  the  third  figure. 
Although  these  regressions  were  not  shown  to  be  significant,  presumably 
because  these  graphs  show  a  larger  scatter  than  had  been  anticipated  from 
previous  studies  when  the  experiment  was  planned,  all  three  graphs  do  seem  to 
show  a  strengthening  o.  anisotropy  during  the  early  stages  of  the  tests  whilst 
the  samples  were  consolidating  with  negligible  lateral  strain,  and  a  weakening 
of  anisotropy  during  the  late  stages  of  the  tests  whilst  the  lateral  strain 
was  increasing. 


5.6  RANDOM  CLUSTERS 

In  addition  to  mapping  domains  defined  by  strong  internal  preferred  orienta¬ 
tion  of  the  clay  plates,  random  clusters  of  clay  plates  were  also  mapped  in 
the  scanning  electron  micrographs.  Fig.  5.8  shows  the  area  of  random  clusters 
vs.  axial  strain.  The  total  random  area  is  fairly  small,  say  5  %  in  comparison 
with  say  95  %  domains.  There  might  have  been  a  trend  concave  upwards,  which 
would  agree  with  the  concave  downward  trend  for  strength  of  anisotropy,  but  if 
so  it  is  small. 


5.7  N'ORMALL  Y-CONSOL I  DATED  DRAINED  TESTS 


SIZES  OF  FEATURES 

Fie.  5.')  .hows  the  average  chor.i  I.-neth  vf  the  areas  mappeu  usine  i'OLMAP 
■'iXial  strain,  the  trend  apparently  be  ins  -oncave  upwards. 

Fig.  5.10  illustrates  the  average  sire  o!  the  iomams  capped  the  scanning 
electron  micrographs  vs.  axial  strain.  Trend  lines  are  shown  lor  horizontal, 
vertical,  and  two  inclined  groups  of  domains.  The  domains  with  horizontal 
preferred  internal  orientation  were  the  largest:  the  vertical  domains  were  the 
•mailest:  tnd  the  inclined  domains  were  intermediate  in  size.  Strain  appears 
to  have  had  negligible  effect  on  the  sizes  of  these  domains;  thus,  both  the 
hypothesis  that  the  domains  would  increase  in  size  and  the  hypothesis  that  the 
domains  would  break  up  and  decrease  in  size  with  shear  appear  to  be 
. 1 iminated. 

Fig.  5.11  shows  the  average  size  of  the  random  clusters  mapped  in  the  scanning 
electron  micrographs  vs.  axial  strain.  These  are  comparable  in  size  to  the 
vertical  domains  and  seem  to  show  a  trend  concave  upwards. 


5.8  POROSITY 

Fig.  5.12  illustrates  the  porosity  of  the  features  mapped  in  the  scanni 
electron  micrographs  vs.  axial  strain.  The  results  here  were  obtained  from  the 
scanning  electron  micrographs  after  Wiener  filtering  and  thresholding  to 
obtain  a  binary  *  sol id-or-pore ’  image.  The  deviations  from  the  trend  lines 
were  relatively  small.  Although  the  trend  line  for  the  total  porosity  obtained 
by  gravimetric  analysis  was  concave  upwards,  see  the  graph  of  volumetric 
strain  above,  the  trend  lines  themselves  were  all  straight.  As  expected,  the 
average  porosity  of  the  random  clusters  was  higher  than  that  of  the  domains. 
Somewhat  unexpectedly,  the  average  porosity  of  the  horizontal  domains  was 
higher  than  that  of  the  vertical  domains;  as  expected,  the  porosities  of  the 
inclined  domains  were  intermediate  between  those  of  the  horizontal  and  verti¬ 
cal  domains.  The  fact  that  all  the  trend  lines  fall  at  the  same  rate  suggests 
that  the  deformation  of  the  samples  is  being  accommodated  within  all  parts  of 
the  structure  at  the  same  rate,  i.e.  that  the  micro-deformations  are  distrib¬ 
uted  uniformly  throughout  the  samples. 


5.9  SUMMARY 

During  the  early  stage  of  shear  deformation,  up  to  8  T,  axial  strain  say:  the 
lateral  strain  was  almost  negligible:  the  preferred  orientation  remained 
horizontal;  the  strength  of  anisotropy  seemed  to  increase;  the  area  of  random 
clusters  might  have  decreased;  the  sizes  of  the  random  clusters  seemed  to 
decrease;  the  sizes  of  the  domains  apparently  remained  constant;  the  sizes  of 
the  areas  mapped  by  optical  microscopy  seemed  to  remain  constant:  and  the 
porosity  of  all  sectors  of  the  structure  decreased  at  the  same  rate.  All  this 
suggests  uniaxial  consolidation  during  the  early  stage  of  shear  deformation, 
with  the  micro-deformations  being  distributed  uniformly  throughout  all  sectors 
of  the  samples. 

During  the  late  siage  of  shear  deformation,  from  8  *  to  16  axial  strain  say: 
the  lateral  strain  increased;  the  preferred  orientation  remained  horizontal. 


5.9  N'ORMTLLY-dQN'SOI . I DATED  DRAINTD . TESTS _ 53 

-it  the  .:r;-r.2th  O!  misccropy  .emeu  ..acre use:  the  uea  c:  .  .nuom  ,:uaU-i'j 
auht  have  :::trtrcn:  •  he  tees  ’  the  random  clusters  cemed  to  uap.-ase:  trie- 
ires  of  the  iomams  apparent  l;.  emuincd  .mutant,  rut  the  -  ires  or  the  areas 
uppeu  by  gt ieai  microscopy  - eemea  'omet imes  to  increase:  md  the  porosity  of 
-ii  ectors  of  the  structure  decreased  at  the  ame  rate.  This  a  train  suggests 
that  the  nncro-de  format  ions  tvere  distributed  uniformly  throughout  ,ai  .ectors 
f  the  samples.  However.  ; t  seems  that  particles  or  domains  ire  rerunning  to 
rotate  away  f  rcm  the  nori/ontai .  out  with  c-quai  quantities  turning  upwards  .md 
downwards :  and  that  these  rotations  caused  'light  iisrupt ion  to  the  structure. 
The  optical  microscopy  - uegests  that  there  may  perhaps  have  been  a  tendency  to 
medium-scale  organization,  with  domains  in  one  part  tending  -ill  to  turn  one 
tay.  and  a ith  domains  in  an  adjacent  part  tending  ill  to  turn  the  other  way. 

in  summary,  therefore,  the  observations  uegest  that  three  different 
mechanisms  were  acting:  one  in  the  ear i y  stage  of  shear  deformation  as  the 
samples  bedded  down:  the  other  in  the  late  stage  of  shear  deformation  as  the 
ainples  began  to  be  disrupted:  and  both  these  mechanisms  were  distinct  from 
the  formation  of  failure  planes. 
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5  10  15 

Axicl  Strain, % 

5.1  Deviator  stress  vs.  Axial  strain 

( ncd ) . 

Axial  Strain,*/* 

2  6  10  U 

18 

Volumetric  strain  vs.  .Axial  strain  (ncd) 


Axiol  Strain,  % 

Fig.  5.4  Orientation  vs.  Axial  strain  (SEM  ncd). 
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Fig. 


Index  of  Anisotropy  vs.  Axial  strain  (SEM  ncd). 


Random  Area.  % 
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5- NORMALLY-CONSOLIDATED  DRAINED  TESTS.  Figure 


Fig.  5.S  Area  of  random  clusters  vs.  Axial  strain  (SEM  ncd). 
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is.  5.10  Size  of  domains  vs.  .Axial  strain  ISEM  ncd). 


Fig.  5.11  Size  of  random  clusters  vs.  Axial  strain  (SEM  ncd) 
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Axial  Strain, •/• 


Fig.  5.12  Porosity  vs.  Axial  strain  (SEM  ncd). 

Fig.  5. 12. a  Horizontal  domains. 

Fig.  5.12.b  H,  I,  V  =  horizontal,  inclined  and  vertical  domains; 

R  =  random  clusters; 


Chapter 


OVER-CONSOLIDATED  DRAINED  TESTS 


6.1  INTRODUCTION 

After  completing  the  normally  consolidated  drained  triaxial  tests,  the  new 
methods  of  micro-analysis  were  then  applied  to  a  series  of  35  over¬ 
consolidated  drained  triaxial  tests  of  Speswhite  Kaolin  after  preparation  in 
the  medium  pressure  filter  presses.  In  brief :  the  slurry  wa^  mixed  at  140  T, 
moisture  content,  i.e.  500  Liquidity  Index:  consolidated  uniaxial  ly  to 
"00-kPa;  consolidated  triaxial  ly  to  ’’00  kPa:  allowed  to  rebound  triaxial  1>  to 
100  kPa:  strained  at  5  %  per  day:  impregnated  with  Ara.dite  AVIS:  and  prepared 
for  microscopy.  All  the  samples  were  examined  in  the  SEM:  but  some  thin  sec¬ 
tions  were  lost  in  the  post,  so  the  optical  microscopy  was  performed  on  IS 
samples.  Apa't  from  rebound,  both  the  normally  consolidated  and  over- 
consolidated  series  of  drained  tests  were  identical. 


6.2  STRESS-STRAIN  BEHAVIOR 

Fig.  6.1  illustrates  the  deviator  stress  vs.  axial  strain  for  a  typical  simple 
from  the  over-consolidated  drained  series.  There  was  a  steady  rise  of  deviator 
stress  to  about  S  °5  axial  strain,  at  about  which  strain  failure  occurred, 
there  being  a  small  amount  of  variability  from  sample  to  sample.  The 
volumeti ic  strain  was  small,  well  less  than  1  %  expansion,  and  there  was  some 
difficulty  in  measuring  it  exactly.  It  should  be  taken  that  it  rose  steadily 
from  the  start,  perhaps  leveling  off  after  about  1-2  %  axial  strain. 


o.3  SCANNING  ELECTRON  MICROSCOPY 

Visual  examination  of  the  scanning  electron  micrographs  showed  that  the  struc¬ 
ture  of  the  samples  was  composes  mainly  of  diffuse  domains,  which  tended  to 
merge  together.  As  indicated  below,  there  were  a  few  random  clusters;  but 
there  were  virtually  no  inter-domain  voids  in  these  samples.  Thus,  the 
hypothesis  that  shear  deformation  is  accommodated  by  slip  of  domains  into  the 
inter-domain  voids  is  irrelevant  here:  but  there  could  still  have  been  slip  of 
the  domains  into  the  looser  and  presumably  weaker  random  areas. 

Failure  planes,  when  they  occurred,  were  slickensided  as  usual  and  cut  through 
the  previously-existing  structure  non-conformably;  however,  attention  here  is 
concentrated  on  the  ’pre-peak’  deformation. 
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o.4  PREFERRED  ORIENTATION 

I'va.  o.2  shows  the  preferred  oriental ion  vs.  ixial  train  obtained  from  the 
polarising  micro-photometric  analysis.  .-MPA.  tie  re  and  in  the  graphs  following, 
each  point  is  the  end  point  of  one  test.  The  results  in  this  figure  are  the 
moduli  of  the  algebraical  average  of  the  five  values  of  orientation  obtained 
from  the  five  fields  measured  by  this  technique,  and  they  are  expressed  as 
inclinations  from  the  horizontal.  Fig.  fa. 3  -hows  the  corresponding  values 
obtained  from  scanning  electron  microscopy.  All  the  values  in  both  graphs  are 
low:  i.e.  the  preferred  orientation  started  close  to  the  horizontal  and 
remained  close  to  the  horizontal  throughout.  Thus,  the  hypothesis  that  the 
particles  would  gradually  rotate  into  the  direction  of  the  failure  plane  is 
e 1 i m i na  t  ed . 


fa. 5  STRENGTH  OF  ANISOTROPY 

Figs.  fa. 4.  o. 5.  and  fa. fa  show  the  strength  of  anisotropy  vs.  axial  strain 
measured  by  PMPA.  POL.MAP,  and  SEM,  respectively.  Each  point  is  the  arithmeti¬ 
cal  average  of  Anisotropy  Index  for  5,  4.  and  24  fields,  respectively. 

In  the  first  and  last  of  these  three  graphs,  the  Anisotropy  Index  starts  at  a 
moderate  value  and  ends  at  a  moderate  value.  This  is  in  accordance  with  the 
hypothesis  that,  during  shear  deformation,  strongly  anisotropic  samples  become 
less  anisotropic  and  weakly  anisotropic  samples  become  more  anisotropic;  these 
samples  appear  to  be  on  the  borderline.  The  values  from  POLMAP  tend  to  be 
higher  throughout. 

All  three  graphs  appear  to  show  a  very  slight  increase  in  anisotropy,  and  the 
standard  deviations  for  SEM  appear  to  show  a  very  slight  decrease  in 
variability  within  the  samples. 


6. fa  RANDOM  CLUSTERS 

Random  clusters  of  clay  plates  were  also  mapped  in  the  scanning  electron 
micrographs.  Fig.  6.7  shows  the  area  of  random  clusters  vs.  axial  strain.  The 
total  random  area  is  fairly  small,  and  there  appears  to  be  a  slight  decrease, 
which  would  agree  with  the  slight  increase  in  anisotropy. 


6.7  SIZES  OF  FEATURES 

Fig.  6.S  shows  the  average  chord  length  of  the  areas  mapped  using  POLMAP  vs. 
axial  strain,  the  trend,  if  any,  being  slightly  upwards  or  concave  down. 

Fig.  6.9  illustrates  the  average  size  of  the  domains  mapped  in  the  scanning 
electron  micrographs  vs.  axial  strain.  Trend  lines  are  shown  for  horizontal 
and  vertical  domains.  There  appears  to  have  been  a  very  slight  increase  in  the 
sizes  of  these  domains. 
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■>.o  POROSITY 

i  in.  o.  10  illustrates  the  porositv  r  -ome  m  the  ieaiures  mapped  .n  the  s  can  - 
tuna  electron  micrographs  vs.  axial  strain.  The  results  nere  were  obtained 
from  the  scanning  electron  micrographs  alter  Viener  I  i  1  terms  and  threshold  ins 
to  obtain  a  binary  ' sol id-or-fxire ‘  tmase.  There  appears  to  have  been  a  very 
slight  decrease  in  porosity  with  strain:  and  this  occurred  more  in  the  •  er t  :  - 
_al  domains  than  in  the  horizontal  ones. 


n.d  SUMMARY 

fp  to  8  ‘T.  axial  strain  say:  the  lateral  strain  was  almost  negligible:  the 
preferred  orientation  remained  horizontal:  the  strength  a I  anisotropy  appeared 
to  increase  slightly:  the  area  of  random  clusters  appeared  to  decrease:  the 
sizes  of  the  domain:  appeared  to  increase  slightly:  and  the  sizes  ol  the  areas 
mapped  by  optical  microscopy  scarcely  varied.  All  this  suggests  very  slight 
bedding  down  of  the  microstructure,  with  groups  of  particles  turning  somewhat 
closer  to  the  horizontal,  and  perhaps  some  oriented  domains  pushing  into  ran¬ 
dom  areas.  The  mechanisms  here  are  very  similar  to  those  in  the  normally  con¬ 
solidated  drained  series  in  the  same  range  of  axial  strain,  and  they  are  again 
•list i net  from  the  formation  of  failure  planes. 
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Axial  Stra>n?  % 


Fig.  6.4  Anisotropy  Index  vs.  Axial  strain  (PMPA  oc ) 


Axial  Strain,  % 

Fig.  6.5  Anisotropy  Index  vs.  Axial  strain  (POLMAP  oc) 


Standard  Deviation  Anisotropy  Index 
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Axial  Strom,  % 


Fig.  6.6 


Index  of  Anisotropy  vs.  Axial  strain  (SEM  oc). 
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Axial  Strain.  % 


Fig.  6 . 7  Area  of  random  clusters  vs.  .\xial  strain  ( SEM  oc ) . 


Axial  Stroin,  % 


Fig.  6.8  Chord  length  of  domains  vs.  Axial  strain  (POLMAP  oc). 
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Axial  Strain,  •/• 


g.  6.9  Size  of  domains  vs.  .^xial  strain  ( SEM  oc ) . 


6  OVER-CONSOLIDATED  DRAINED  TESTS.  Figures 
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Axial  Strain,  % 


Fig.  6.10  Porosity  vs.  .Axial  strain  ( SEM  oc ) . 
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NORMALLY-CONSOLIDATED  L'NDRA  I  N  E  D  TESTS 


' . i  INTRODUCTION 

The  third  aeries  of  tests  to  he  considered  is  a  series  of  35  normally  con¬ 
solidated  undrained  triaxial  tests  of  Speswhite  kaolin  after  preparation  in 
the  original  low  pressure  filter  presses.  The  kaolin,  although  nominally 
identical  to  that  used  in  the  other  two  series,  was  slightly  different.  In 
brief:  the  slurry  was  mixed  at  120  moisture  content,  i.e.  2.56  %  Liquidity 
Index:  consolidated  uniaxial ly  to  100  kPa:  consolidated  triaxial ly  to  650  kPa: 
itrained  at  0.15  mm/min  under  50  kPa  backpressure:  impregnated  with  Yestopal 
'.V:  and  prepared  for  microscopy. 


".2  STRESS-STRAIN  BEHAVIOR 

fig.  ".1  illustrates  the  deviator  stress  vs.  axial  strain  for  three  typical 
samples  from  the  normally  consolidated  undrained  series,  the  results  being 
much  as  expected,  i.e.  a  steady  rise  of  deviator  stress  to  a  peak  value  at  say 

S  T.  axial  strain,  followed  by  a  more-or-less  steady  stress  to  say  16  %  axial 

strain,  at  about  which  strain  failure  occurred.  Confidence  Limits  for  the 

whole  series  are  shown  on  the  diagram:  the  spread  of  these  at  high  strains  is 

due  in  part  to  the  smaller  number  of  samples  which  were  tested  to  these  high 
strains. 


T3  VISUAL  EXAMINATION 

Visual  examination  in  the  polarizing  microscope  showed  that  a  certain 
anisotropy  had  developed  during  the  preparation  of  the  filter  cake  and  that 
most  of  the  particles  had  turned  into  the  horizontal  direction:  but,  after 
triaxial  consolidation,  a  blotchy  appearance  could  be  seen  between  crossed 
polars,  and  the  particles  were  arranged  somewhat  randomly. 

Visual  examination  of  the  scanning  electron  micrographs  showed  that  the  struc¬ 
ture  of  the  samples  was  composed  mainly  of  diffuse  domains,  which  tended  to 
merge  together.  As  indicated  be  low, there  were  a  few  random  clusters:  but  there 
were  virtually  no  inter-domain  voids  in  these  samples.  Thus,  the  hypothesis 
that  shear  deformation  is  accommodated  by  slip  of  domains  into  the  inter- 
domain  voids  is  irrelevant  here:  but  there  could  still  have  been  slip  of  the 
domains  into  the  looser  and  presumably  weaker  random  areas. 

Failure  planes,  when  they  occurred,  were  slickensided  as  usual  and  cut  through 
the  previously-existing  structure  non-conformably:  however,  attention  here  is 
concentrated  on  the  ’pre-peak’  deformation. 
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The  measurement •'  u is cussed  below  were  t .tk.cn  -on  fields  >a  f.'i.ows: 

I'MPA  !.5  mm  diameter 

POLMAP  0.4"  mm  square 

SEM  0.025  mm  square 


~ . 4  PREFERRED  ORIENTATION 

Fig.  '.2  shows  the  preferred  orientation  vs.  axial  strain  obtained  from  the 
polarizing  nucro-photomet r ic  analysis.  PMPA.  Here  and  in  the  graphs  following, 
each  point  is  the  end  point  of  one  test.  Here  and  in  the  following  graphs, 
quadratic  regressions  have  been  added.  Although  these  regressions  were  not 
-.hown  to  be  significant,  presumable  because  these  graphs  show  a  larger  scatter 
than  had  been  anticipated  from  previous  studies  when  the  experiment  was 
planned,  many  of  the  graphs  show  the  same  trends.  The  results  in  Fig.  ~.2  are 
the  moduli  of  the  algebraical  average  of  the  five  values  of  orientation 
obtained  from  the  five  fields  measured  by  this  technique,  and  they  are 
expressed  as  inclinations  from  the  horizontal.  The  values  for  the  averages  of 
the  five  moduli  were  slightly  higher,  and  their  trend  was  slightly  upwards, 
rig.  ~.3.  The  trends  for  the  sets  of  four  measurements  obtained  from  POLMAP 
are  shown  in  Fig.  ~.j.  Fig.  ".4  shows  the  values  obtained  from  scanning  elec¬ 
tron  microscopy  using  intensity  gradient  analysis  followed  by  consistency 
ratio  smoothing.  PICCON;  and  Fig.  ~.5  shows  the  values  obtained  from  the 
porosity/density  maps  by  using  intensity  gradient  analysis  on  these  maps.  PIC- 
POR.  All  the  values  in  these  four  graphs  are  low;  i.e.  the  preferred  orienta- 
t ion  started  close  to  the  horizontal  and  remained  close  to  the  horizontal 
throughout.  Thus,  the  hypothesis  that  the  particles  would  gradually  rotate 
into  the  direction  of  the  failure  plane  is  eliminated;  on  the  contrary,  the 
trend  being  concave  up  suggests  that  during  the  early  stage  of  the  tests,  up 
to  3  T,  axial  strain  say,  the  particles  are  turning  towards  the  horizontal,  and 
that  during  the  late  stage  of  the  tests,  from  S  -  16  T,  axial  strain  say.  some 
of  the  particles  may  be  turning  away  from  the  horizontal. 


~ . 5  STRENGTH  OF  ANISOTROPY 

Figs.  ^.6,  ~.S.  and  7.9  show'  the  strength  of  anisotropy  vs.  axial  strain 

measured  by  PMPA,  POLMAP.  PICCON  and  PICPOR.  respectively.  Each  point  in  the 
first  two  graphs  is  the  arithmetical  average  of  .Anisotropy  Index  for  5  and  4 
fields,  respectively;  and  each  point  in  the  other  two  graphs  is  the  composite 
value  of  Consistency  Ratio  for  25  fields. 

In  all  four  graphs,  the  Anisotropy  Index  or  Consistency  Ratio  starts  at  a  mod¬ 
erate  value  and  ends  at  a  moderate  value.  This  is  in  accordance  with  the 
hypothesis  that,  during  shear  deformation,  strongly  anisotropic  samples  become 
less  anisotropic  and  weakly  anisotropic  samples  become  more  anisotropic;  these 
samples  appear  to  be  on  the  borderline. 

All  four  graphs  do  seem  to  show  a  strengthening  of  anisotropy  during  the  early 
stages  of  the  tests,  and  a  weakening  of  anisotropy  during  the  late  stages  of 
the  tests. 

The  values  from  POLMAP  were  interesting  in  that  the  between-samples 
variability  seemed  to  increase  with  strain;  this  point  is  discussed  further 
below. 
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.\nother  interest ina  feature  :  that  the  -trenath  of  mi  sot  ropy  based  on  the 
density  maps  of  PICFOR  is  mucn  lower  than  that  rased  on  the  t lent at  ion  maps 
of  PICCOX:  this  raises  the  suggestion  that  the  m isotropic  mechanical  proper¬ 
ties  of  the  samples  may  be  due  to  the  anisotropy  of  density  rattier  than  to  the 
anisotropy  of  orientation.  This  point  deserves  further  study. 


'.o  SIZES  OF  FEATURES 

Tig.  '.10  shows  the  average  chord  length  of  the  areas  mapped  using  POLMAP  vs. 
axial  strain,  the  trend  apparently  being  concave  downwards:  however,  the 
observers  were  not  altogether  convinced  that  the  strong  downward  trend  of  the 
quadratic  regression  did  represent  what  they  had  seen. 

Figs.  '.11  and  '.12  illustrate  the  average  size  of  the  features  mapped  in  the 
'.canning  electron  micrographs  by  PICCON  and  PICPOR.  respectively.  Strain 
appears  to  have  had  negligible  effect  on  the  sizes  of  these  features:  thus, 
both  the  hypothesis  that  the  domains  would  increase  in  size  and  the  hypothesis 
that  the  domains  would  break  up  and  decrease  in  size  with  shear  appear  to  be 
el i m i na  t  ed . 


'.'  POROSITY 

The  relationship  between  porosity  and  strength  of  anisotropy  was  investigated 
on  a  pixel-by-pixel  basis.  For  all  the  pixels  in  all  25  images  of  each  set. 
the  local  values  of  Consistency  Ratio  and  of  porosity  were  taken  from  the 
smoothing  filter:  and  regressions  of  one  on  the  other  were  calculated.  It  was 
found  that  the  more  strongly  oriented  domains  were  the  denser. 

A  trial  analysis  by  hand  on  a  few  micrographs  suggested  that  the  denser  parts 
of  the  samples  contained  larger  domains  than  did  the  looser  parts. 

Taken  together,  these  observations  suggested  the  possibility  that  the  micro¬ 
deformations  were  concentrated  in  the  looser  parts  of  the  samples,  the  denser 
parts  remaining  unaltered. 


7.3  BETWEEN- SAMPLES  VARIABILITY 

Figs.  '.13  and  '.14  show  the  between-samples  standard  deviations  for  strength 
of  anisotropy  and  for  chord  sizes  obtained  from  the  optical  microscopy  calcu¬ 
lated  by  placing  the  samples  in  groups  with  approximately  equal  strains  cor¬ 
responding  to  the  vertices  on  the  graphs.  These  graphs  suggest  that  the 
variability  between  samples  increased  as  strain  increased. 


7.9  SUMMARY 

During  the  early  stage  of  shear  deformation,  up  to  S  %  axial  strain  say:  the 
preferred  orientation  remained  close  to  the  horizontal:  the  strength  of 
anisotropy  seemed  to  increase:  the  sizes  of  the  domains  apparently  remained 
constant;  and  the  sizes  of  the  areas  mapped  by  optical  microscopy  sometimes 
increased. 
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'•ur-:;u  t  he  ate  taee  of  >hear  de  format ;  or; .  :  ron;  to  '  axiui  .train  -u\: 
t  he  preferred  •.  ri  dilation  remained  horiconta. .  out  t  treneth  ii  in isotropy 
- eemed  to  decrease:  the  sizes  at"  the  domains  apparent  iy  remained  constant  .  out 
'.lie  i zes  .:f  the  ireas  mapped  by  optical  microscopy  -cemed  sometimes  to 
increase:  and  the  be  tween-samples  variability  apparent !y  i ncreaseo . 

riii.s  suggests  that  particles  or  domains  are  '-e  ginning  to  rotate  .‘.way  from  the 
horizontal,  but  with  equal  quantities  turning  upwards  and  downwards :  and  that 
these  rotations  caused  slight  disruption  to  the  structure.  The  optical  micro¬ 
scopy  suggests  that  there  may  perhaps  have  been  a  tendency  to  medium-scale 
organization,  with  domain-,  in  one  part  tending  all  to  turn  one  way.  md  -with 
domains  in  an  adjacent  part  tending  all  to  turn  the  other  way. 

In  summary,  therefore,  the  observations  suggest  that  three  different 
mechanisms  were  acting:  one  in  the  early  stage  of  shear  deformation  as  the 
samples  bedded  down:  the  other  in  the  late  stage  of  shear  deformation  as  the 
samples  began  to  he  disrupted:  and  both  these  mechanisms  were  distinct  from 
the  formation  of  failure  planes.  There  may  also  have  been  some  form  of  struc¬ 
tural  instability,  such  that  samples  which  were  nominally  identical  actually 
behaved  somewhat  differently. 

I'he  questions  of  the  relative  importance  of  anisotropy  of  density  and 
anisotropy  of  orientation  and  of  the  relative  roles  of  the  denser  and  looser 
parts  of  the  structure  might  well  form  the  subject  of  further  tudies. 
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Fig.  '.2  Orientation  vs.  Axial  strain  (PMPA  ncu ! 
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Fig.  ~.3  Orientation  vs.  .Axial  strain  (PMPA  and  POLMAP 
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Fig.  ~.b  .Anisotropy  Index  vs.  .Axial  strain  iPMPA  neu ) . 
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■’.?  .Anisotropy  Index  vs.  Axial  strain  (POLMAP  ncu). 
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Fig.  7.10  Size  of  patches  vs.  Axiai  strain  (POLMAP4  ncu) 
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Fig.  ".II  Size  of  domains  vs.  Axial  strain  ! SEM  PICCON  ncu ) 
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Fig.  7.12  Size  of  features  vs.  Axial  strain  (SEM  PICPOR  ncu) 
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Fig.  ".13  Be tween-samples  standard  deviation  of  anisotropy  vs.  Axial  strain 

( PMPA .  POLMAP.  POLMAP4  ncu  ) . 


Fig.  7.14  Between-samples  standard  deviation  of  size  of  patches  vs.  Axial 

strain  (P0LMAP4  ncu). 
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